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ing  cylinder  with  a  thin  dielectric/ferrite  coating.  The  formal  integral  representation  for  the  field  is  evaluated  asymptotically 
via  the  stationary  phase  method  for  the  field  in  lit  region  and  via  the  residue  theorem  for  the  field  in  shadow  region.  The 
obtained  solution  is  cast  in  the  format  of  the  geometrical  theory  of  diffraction  (GTD).  In  the  lit  region,  the  geometrical  optics 
(GO)  solution  consists  of  the  direct  incident  ray  and  the  reflected  ray.  The  ray  picture  associated  with  the  reflected  ray 
for  the  coated  cylinder  remains  the  same  as  that  for  the  conducting  cylinder  without  coating.  The  reflection  coefficient  for 
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the  uniform  solution.  Numerical  data  for  the  essential  transition  integral  is  deduced,  via  a  heuristic  approach,  from  the 
exact  eigenfunction  solution  together  with  the  GTD  solutions  for  the  coated  cylinder.  The  solution  is  convenient  for  the 
engineering  applications  due  to  its  simple  ray  format.  Numerical  results  obtained  from  the  UTD  solution  for  the  cylinder 
with  thin  coating  show  excellent  agreement  with  the  exact  eigenfunction  results. 
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CHAPTER  I 


INTRODUCTION 


1.1  Outline 

The  purpose  of  this  research  is  to  develop  an  approximate  asymptotic  high- 
frequency  solution  for  the  electromagnetic  scattering  from  a  convex  conducting 
surface  with  dielectric/ferrite  coating.  For  the  bare  perfectly  conducting  case,  an 
efficient  uniform  GTD  solution  has  been  developed  and  used  in  many  engineer¬ 
ing  apphcations.  However,  no  such  'iolution  exists  for  the  conductor  coated  with 
dielectric/ferrite  material. 

This  research  is  motivated  by  the  need  for  understanding  the  effects  of  loading 
on  the  electromagnetic  scattering.  In  some  engineering  problems  there  is  a  need 
to  control  the  scattered  fields  from  conducting  objects  such  as  the  reduction  of  the 
radar  cross  section  (RCS)  in  the  design  of  air- borne  vehicles.  An  efficient  method  to 
control  the  EM  scattering  can  be  achieved  by  coating  of  a  dielectric/ferrite  material 
on  the  conducting  scatterers  without  having  to  deform  the  physical  structure  of 
the  scatterer.  The  research  in  this  report  is  focused  on  the  development  of  an 
effficient  analytical  tool  to  calulate  the  scattered  field  from  a  coated  surface  and 
the  effects  of  various  coating  materials  on  the  scattered  field. 

It  is  desired  to  ultimately  develop  the  uniform  geometrical  theory  of  diffrac¬ 
tion  (UTD)  for  an  general  convex  surface  with  material  coating.  As  a  first  step,  we 
examine  the  field  scattered  by  an  infinite  conducting  circular  cylinder  with  dielec¬ 
tric/ferrite  coating  when  it  is  illuminated  by  a  parallel  infinite  electric  or  magnetic 
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line  source.  Since  both  the  line  source  and  the  cylinder  are  parallel  and  have 
infinite  lengths,  the  scattering  from  the  cylinder  is  essentially  a  two-dimensional 
problem.  Once  the  solution  for  this  canonical  problem  is  obtained,  it  can  be  easily 
extended  to  calculate  the  field  scattered  by  a  general  convex  coated  surface.  The 
scattering  from  the  cylinder  with  surface  impedance  is  also  treated  in  this  research. 
This  problem  is  of  interest  in  that  it  provide  the  solution  of  scattering  field  from 
a  cylinder  with  finite  conductivity. 

The  organization  of  this  report  is  as  follows.  Chapter  II  gives  the  angular 
eigenfunction  series  solutions  for  the  impedance  and  the  coated  cylinder.  It  is  noted 
that  the  eigenfunction  solution  is  slowly  convergent  when  the  radius  of  cylinder  is 
large  in  terms  of  the  wave  length.  Even  though  the  eigenfunction  solution  is  not 
the  goal  of  this  research,  it  is  important  because  it  provides  a  means  of  checking 
the  accuracy  of  the  ray  solution  which  is  developed  in  later  Chapters. 

In  Chapter  III,  high  frequency  asymptotic  solutions  are  obtained  from  the 
asymptotic  evaluation  of  the  field  expression  in  the  integral  form.  The  solutions 
obtained  are  cast  in  the  format  of  the  geometrical  theory  of  differation  (GTD) 
as  suggested  by  Keller  [1,2].  In  the  lit  region,  geometrical  optics  fields  (i.e.  the 
incident  and  the  reflected  field  ^  are  obtained  via  the  stationary  phsise  method. 
In  the  shadow  region,  the  diffracted  field  is  obtained  from  the  evaluation  of  the 
integral  via  Cauchy’s  residue  theorem.  It  is  noted  that  the  ray  picture  for  the 
impedance  and  the  coated  cylinder  remains  the  same  as  that  for  the  perfectly 
conducting  cylinder.  Thus  the  Keller’s  GTD  ray  format  given  for  the  conducting 
cylinder  can  be  used  with  necessary  modifications  due  to  the  boundary  conditions 
on  the  cylinder  surface.  Even  though  the  solution  obtained  is  accurate  in  the  deep 
lit  region  and  in  the  deep  shadow  region,  it  fails  in  the  shadow  boundary  transition 
region.  The  reasons  for  this  failure  are  also  analyzed  in  this  Chapter. 

Chapter  IV  examines  some  characteristics  of  the  creeping  waves  which  prop- 
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agate  azimuthally  on  the  cylinder  surface.  The  propagation  constants  and  the 
attenuation  constants  of  the  creeping  wave  modes  are  determined  from  the  zeros 
of  the  transcendantal  equation  which  characterizes  the  impedance  or  the  coated 
cylinder.  For  the  conducting  cylinder,  the  contribution  from  the  first  creeping 
wave  mode  to  the  diffraction  is  always  dominant  over  other  higher  modes.  How¬ 
ever  this  is  not  always  true  for  the  impedance  and  the  coated  cylinder.  Change 
of  the  dominance  of  the  creeping  wave  mode  according  to  the  change  of  the  rele¬ 
vant  parameters  is  illustrated  by  the  some  examples.  The  behavior  of  the  creeping 
waves  for  a  very  large  cylinder  is  compared  with  the  surface  waves  of  the  grounded 
planar  slab  with  the  surface  impedance  or  the  dielectric/ferrite  coating. 

In  chaper  V,  the  failure  of  the  GTD  solution  in  the  shadow  boundary  tran¬ 
sition  region  is  corrected  by  the  heuristic  extension  of  the  UTD  solution  for  the 
conducting  cylinder  to  the  impedance  and  the  coated  cylinder.  For  the  impedance 
and  the  coated  cylinder,  the  transition  function  P{T,q),  which  is  well  tabulated 
for  the  conducting  case,  is  not  available.  Instead  of  direct  numericad  evaluation  of 
the  integreil,  a  heuristic  method  is  used  to  obtain  the  numerical  data  for  the  tran¬ 
sition  function  P{x,q).  The  universal  feature  of  this  transition  function  is  tested 
for  some  examples.  Using  the  transition  function  obtained,  the  field  patterns  are 
calculated  and  the  results  are  compared  with  exact  eigenfunction  solutions. 

Chapter  VI  contains  the  conclusion  for  this  research. 

1.2  Literature  survey 

For  the  perfectly  condvs  ing  case,  the  asymptotic  evaluation  of  the  field  scat¬ 
tered  by  a  cylinder  or  a  general  convex  surface  has  received  much  attention  and 
been  studied  extensively.  Such  solutions  are  usually  cast  in  the  format  of  the  ge¬ 
ometrical  theory  of  differation  (GTD)  as  suggested  by  Keller  [1,2].  One  of  more 
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t-:.  £nt  works  on  this  problem  is  that  of  Pathak  [3].  He  made  a  significant  contri¬ 
bution  to  this  area  by  developing  an  approximate  uniform  geometrical  theory  of 
diffraction  (UTD)  for  a  perfectly  conducting  cylinder  and  arbitrsu-y  convex  surface 
with  great  success. 

The  problem  of  scattering  by  the  impedance  or  the  coated  cylinder  has  also 
received  attention  for  decades.  Scattered  field  by  the  impedance  cylinder  has  been 
studied  by  Wait  [17],  Streifer  [18]  and  Logan  [14].  The  main  concern  of  their  work 
was  in  ground  wave  propagation  excited  a  dipole  on  the  earth  which  has  been 
modeled  by  a  cylinder  with  surface  impedance.  In  their  work,  the  propagation 
constants  for  the  creeping  wave  of  the  impedance  cylinder  has  been  obtained  with 
the  assumption  that  the  Hankel  functions  in  the  formal  solution  may  be  approxi¬ 
mated  by  their  Fock  type  Airy  functions.  It  is  noted  that  this  assumption  is  not 
valid  for  arbitray  orders  and  arguments  of  the  Hankel  functions.  Recently,  Wang  [4] 
obtained  more  accurate  numerical  data  for  the  propagation  constant  of  the  creep¬ 
ing  wave,  the  Regge  poles,  and  the  natural  frequencies  of  an  impedance  cylinder 
by  using  both  the  Debye  and  the  Watson’s  approximation  of  Hankel  functions. 
The  creeping  wave  interpretation  for  the  resonance  and  the  correlations  between 
the  resonance,  the  Regge  poles,  and  the  natural  frequencies  are  also  discussed. 
Most  recent  work  on  the  scattering  from  the  impedance  cylinder  has  been  done 
by  Wang  and  Kim  [5]  who  has  developed  uniform  GTD  solution  for  this  problem. 
The  accuracy  of  his  results  has  been  confirmed  by  excellent  agreement  with  exact 
eigenfunction  solutions. 

For  the  coated  cylinder,  Tang  (1957)  has  found  an  angular  eigenfunction  series 
solution  for  the  back  scattered  field  and  has  presented  numerical  and  experimental 
results  when  the  incident  wave  length  A©  is  comparable  with  radius  of  the  cylin¬ 
der.  The  propagation  constant  for  the  first  creeping  wave  mode  on  the  surface  of 
cylinder  with  lossless  dielctric  coating  was  first  obtained  by  Elliott  [19]  by  using 
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some  asymptotic  representations  of  the  various  necessary  Bessel  and  Hankel  func¬ 
tion.  Based  upon  Poisson’s  sum  formula,  Kodis  |21]  has  obtained  an  asymptotic 
solution  for  the  scattered  field  in  terms  of  the  difference  between  the  scattered  field 
of  the  coated  cylinder  and  of  a  perfectly  conducting  cylinder.  In  his  solution,  the 
scattered  field  is  expressed  in  term  of  several  different  rays  by  tracing  the  compli¬ 
cated  internal  process  of  the  field  propagation  in  the  dielectric  coating.  Because 
of  difficulty  in  tracing  the  optical  rays  inside  the  coating  material,  his  solution  is 
too  complicated  for  numerical  calculations.  Particularly  interesting  is  the  work  by 
Helstrom  [22],  who  has  obtained  a  ray  optical  solution  by  ausymptotic  evaluation 
of  the  formal  solution.  A  trajectory  of  the  propagation  constants  of  the  creeping 
wave  has  been  sketched.  However,  he  failed  to  give  exact  numerical  data. 

By  employing  the  same  optical  rays  as  given  by  Kodis,  Rao  and  Hamid  [23] 
invesigated  GTD  solution  for  this  problem.  No  numerical  results  were  presented 
in  their  paper.  The  backscattering  width  of  a  dielectric  coated  cylinder  hew  been 
obtained  via  a  high  frequency  ray  solution  by  the  recent  work  of  Wang  [6].  Most 
recently,  Paknys  and  Wang  [7,8]  have  developed  an  approximate  solution  for  the 
surface  fields  induced  on  a  coated  cylinder  by  a  magnetic  or  electric  line  source  on 
the  same  surface.  Even  though  the  final  goal  of  this  work  is  different  from  theirs, 
some  of  their  results  associated  the  propagation  constants  of  the  creeping  waves 
are  utilized  in  this  work. 
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CHAPTER  II 


EIGENFUNCTION  SERIES  SOLUTION 


In  this  Chapter,  the  eigenfunction  series  solutions  for  the  field  scattered  by 
the  infinite  circular  conducting  cylinder  with  an  impedance  boundary  surface  or  a 
dielectric/ferrite  coating  will  be  presented.  In  the  eigenfunction  solution,  the  field 
expression  is  given  as  the  sum  of  each  angular  eigenfuntion  mode  which  satisfies 
the  boundary  conditions  of  a  given  canonical  problem.  The  angular  expansion  of 
eigenfunction  modes  is  an  appropriate  solution  when  the  radius  of  the  cylinder  is 
small  in  terms  of  the  wavelength  where  the  asymptotic  solution  cannot  b*  applied. 
However,  for  the  cylinder  with  large  radius  it  is  impractical  to  compute  the  field 
from  the  eigenfuntion  expansion  because  it  is  very  slowly  convergent. 

Even  though  the  eigenfunction  solution  is  not  the  final  goal  of  this  research,  it  is 
important  because  the  formulation  of  a  high  frequency  ray  solution  begins  with  the 
field  expression  in  the  eigenfunction  series  form.  Since  the  eigenfunction  solution 
is  based  on  exact  boundary  conditions,  it  is  considered  to  be  the  exact  solution  of 
the  given  problem.  Hence,  it  can  be  also  used  as  a  means  to  check  the  validity 
of  the  high  frequency  ray  solution  which  will  be  developed  in  later  Chapters.  In 
addition,  the  eigenfunction  solution  provides  the  basis  to  obtain  the  numerical  data 
for  the  transition  function  Pii^q)  associated  with  the  uniform  GTD  solution  in 
the  vicinity  of  the  shadow  boundary. 
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2.1  The  canonical  problem 


The  geometry  of  the  problem  for  an  impedance  and  a  coated  cylinder  is  il¬ 
lustrated  in  Figure  2.1  (a)  and  (b)  respectively.  The  impedance  cylinder  or  the 
coated  cylinder  is  illuminated  by  a  parallel  line  source  located  at  4>‘).  The 

impedance  cylinder  is  an  impenetrable  cylinder  with  an  impedance  boundary  sur¬ 
face  and  the  coated  cylinder  has  a  material  coating  on  the  conducting  core  cylinder. 
The  coating  has  a  thickness  d  =  b  — a  and  is  homogeneous  with  permittivity  cj  and 
permeability  /ij  ,  both  of  which  may  be  complex.  The  source  is  either  an  electric 
current  filament  /  (TMz  case)  or  a  magnetic  current  filament  K  {TEg  case).  The 
solution  to  be  investigated  is  the  total  field  at  P(/>,  ^)  from  the  line  source  which 
radiates  in  the  presence  of  a  parallel  cylinder  with  an  impedance  boundary  surface 
or  with  a  material  coating.  Both  the  source  and  the  field  point  are  assumed  to  be 
located  outside  cylinder.  Since  both  the  line  source  and  the  cylinder  are  assumed 
to  have  infinite  lengths  and  to  be  parallel  to  the  z-axis,  the  scattering  from  the 
cylinder  is  essentially  a  two-dimensional  problem.  The  two-dimensional  line-source 
Green’s  function  G{p,p')  satisfies  the  equation 


(  V?  +  fco  )  0(p,?)  =  -8{p-  P)  (2.1) 


in  the  domain  b  <  {p,p')  <  <x>  \  <  27r  where  is  two-dimensional 

Laplacian  operator.  Thus  Equation  (2.1)  can  be  rewritten  as 


lA  A  i.  1.2 

p  dp  ^  dp^  p^  d(f>^  ^  ° 


G{p,P) 


~6(p  -  p’)  8{4>  -  4!) 
P 


(2.2) 


Without  loss  of  generality,  we  can  set  =  0  and  0  <  (^  <  rr  for  convenience.  Once 
the  solution  G{p,p')  for  Equation  (2.2)  subject  to  the  boundary  conditions  on 
the  cylinder  surfaces  and  the  Sommerfeld  radiation  condition  as  p  — »  00  is  found, 
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the  field  components  at  P{p^<f>)  can  be  readily  obtained.  For  excitation  by  a  line 

mm  w* 

source  of  electric  current  I  {TMz  case),  the  field  components  of  E  and  Ji  outside 
the  cylinder  are  given  as  (for  time  variation) 


Ez  =  Hz  —  0 

Ep-0, 

=  0, 

For  the  magnetic  line  source  case  (TEz  case),  the  field  components  are  given  as 


Hz  =  -ju}€oKG,  Ez  =  0 

H,  =  a,  =  (2.4) 

=  0,  E^  =  K  ^ 


y 


(a) 


Figure  2.1:  The  geometry  of  the  problem 
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2.2  Impedance  Cylinder 


In  this  section,  the  eigenfunction  solution  for  Equation  (2.2)  subject  to  the 
boundary  condition  of  the  impedance  cylinder  will  be  obtained.  On  the  sur- 
face  of  the  impenetrable  impedance  cylinder,  the  field  G(p,p')  shovild  satisfy  the 
impedance  boundary  condition  given  as 

=  0  (2.5) 

op 

where  C  represents  either  normalized  surface  impedance  (TE,  case)  or  normalized 
admittance  (TMg  case).  G  is  assumed  to  be  constant  around  the  surface  of  the 
cylinder. 

C  =  Cs  =  ^  for  TM^  :  C  =  C*  =  ~  for  TE,  (2.6) 

Zj  Zo 

where  Zo  is  the  intrinsic  impedance  of  free  space  (i.e.  Zo  =  )  and  Z$  is  the 

constant  impedance  on  the  cylinder  surface.  The  subscript  s  and  h  in  C  denote 
TMg  (soft)  and  TEx  (hard)  respectively. 

2.2.1  TMx  case 

At  the  field  point  P{p,  <f>),  the  total  field  consists  of  the  two  terms;  the  incident 
field  Go{p,4>)  and  the  scattered  field  G,{p^<^). 

G{p,P)  =  GM)  =  Go{p,ii>)  +  G.{p,<l>)  (2.7) 

When  the  line  source  is  an  electric  current  filament  /,  the  incident  electric  field 
is  parallel  to  the  cylinder  axis  which  is  taken  as  the  z-axis.  From  the  free  space 
Green's  function  for  two-dimensional  line  source,  the  z-polarized  incident  electric 
field  Go{p,  <i>)  is  given  as 
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(2.8) 


*  •  f2}  — •  ** 

Using  the  additioneJ  therem  for  Ho  '{ko\p  —  p*\),  Go(p,^)  can  be  expanded  in  a 
series  form  given  by 

Go(p,4')  =  ^  E  d%oP>)MKp<) (2.9) 

where  p>  =  p'  ,  P<  ~  P  for  P  <  p\  P>  =  P  •>  P<  —  p'  for  p  >  p'  .  It 
is  pointed  out  that  the  scattered  field  from  an  infinite  parallel  cylinder 

by  z-polarized  incident  field  Go{p,<f>)  is  also  z-polarized.  Considering  the  Equation 
(2.2)  and  the  radiation  condition  at  infinity,  we  assume  that  the  scattered  field 
from  the  cylinder  can  be  expressed  as 

G,ip,<i>)  =  ^  E  A,{n)  Hl?\kop<)Hi?\koP>)  (2.10) 

The  total  field  at  P{p,4>)  is  the  summation  of  the  incident  field  and  the  scattered 
field.  Thus,  G{p,  4>)  is  given  as 

G{p,<i>)  =  ^  E  H^n\koP>)\Jn{koP<)^A,{n)H^n\koP<)]e~^^^  (2.11) 


The  unknown  scattering  coefficient  A5(n)  in  Equation  (2.11)  can  be  easily  deter¬ 
mined  from  the  boundary  condition  given  in  Equation  (2.5).  Thus,  A,(n)  is  given 
as 


Jn{koh)~jCMkoh) 

H^n^\koh)  -  jC>H^n\koh) 


(2.12) 


The  primes  in  Equation  (2.12)  represent  derivatives  with  respect  to  the  arguments. 
For  the  unification  of  notation,  all  Bessel  functions  included  are  replaced  by  Han- 
kel  functions  Hn\x)  and  Hj^^{x).  From  the  well-known  relationships  of  Bessel 
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functions  and  Hankel  functions,  Equation  (2.11)  can  be  rewritten  as 


=  £  Hl?hkoP>)[HL'hkoP<)-i-F.in)Hi^HkoP<)]  i2.U) 


where 


Bi^^\kob)-jCsHi^kkob) 


(2.14) 


A  plane  wave  incident  is  a  limiting  case  when  the  source  is  located  at  infinity  (i.e. 
p'  —*  oo).  The  total  field  for  plane  wave  incidence  can  be  readily  obtained  by 
replacing  Hn\kop')  by  its  large  argument  form.  Thus,  G{p,^)  for  the  plane  wave 
incidence  can  be  expressed  as 


<?(p,<^)-^  Y,  frfW)  +  i".(n)rf^(*c,p)l  (2.15) 

^  n=-oo  '•  * 

where  Vq  is  the  incident  field  at  reference  point  which  is  chosen  at  center  of  cylinder. 
2.2.2  TEz  case 

When  the  line  source  is  a  magnetic  current  filament  /f,  the  source  radiates  a 
z-polarized  incident  magnetic  field  Cto(p,^)  which  is  also  given  as 


Oo(pA)  =  -ja?\K\p-?\) 


(2.16) 


Using  a  similar  procedure  as  for  the  TMz  case,  the  total  field  at  P{p,<^)  is  obtained 
as 


0(pA)  =  -[  E  +  (2.17) 

^  n=-oo  ^  •* 
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The  scattering  coefficient  Af^{n)  subject  to  the  boundary  condition  in  Equation 
(2.5)  is  given  as 

Jn(kob)-jChJnikob) 


Ah{n)  =  - 


Hi?^\kob)  -  jChH[^\kob) 


(2.18) 


If  we  replace  Jn(*)  and  Nn{^)  with  J{n\x)  and  Hn^\x),  Equation  (2.17)  can  be 
rewritten  as 

G{p,4>)^Y  £  4^W>)[rf^(A:aP<)  +  FA(n)rf)(W<)]e--”‘^  (2.19) 


r{2). 


where 


Fhin)  = 


ffj^hkeb)  -  jChHi^\kob) 
H^nhkob)-jChH^n\kob) 


(2.20) 


For  the  case  of  plane  wave  incidence,  the  total  field  (?(p,  <j>)  is  given  as 


^  n=-oo  *• 


(2.21) 
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2.3  Coated  Cylinder 


Following  the  same  procedure  as  for  the  impedance  cylinder,  the  eigenfunc¬ 
tion  solution  for  the  field  scattered  by  the  infinity  conducting  cylinder  with  dielec¬ 
tric/ferrite  coating  for  both  TE^  and  TMz  cases  is  obtained  in  this  section.  For  the 
coated  cylinder,  it  should  be  noted  that  fields  can  penetrate  into  the  air /coating 
interface  (p  =  b)  on  the  given  boundary  condition  that  tangential  components  of 
both  the  electric  and  magnetic  fields  should  be  continuous  at  the  interface. 

At  the  surface  of  the  core  conducting  cylinder,  fields  should  satisfy  the  soft  or  hard 
boundary  conditions  according  to  the  polarization  of  the  incident  fields. 

G{a,4>)  =  0  ;for  soft  boundary  {TMt  case)  (2.22) 


=  0  ;/or  hard  boundary  {TEz  case)  (2.23) 

p=a 

The  above  boundary  conditions  simply  mean  that  on  the  surface  of  the  conducting 
core  cylinder,  the  tangential  component  of  the  electric  field  vanishes  for  both  the 
TMz  and  the  TEz  cases. 

t 

2.3.1  TMz  case 

Since  the  incident  field  from  the  source  is  independent  of  the  scattercr,  the 
expression  for  the  incident  field  Go{p,4>)  remains  unchanged. 

.  oo 

=  ^  E  Hl?\koP>)MkoP<)e-j^^  (2.24) 

n=— oo 

The  scattered  field  Gt{p,4>)  outside  the  cylinder  (i.e.  p  >  b)  can  be  expressed  as 


dG{p,<i,) 

dp 


14 


(2.25) 


G,{p,4>)^~  E  Mn)H^n\koP>)H^n\koP<)e'^^^ 

^  n=— oo 


Thus  the  total  field  outside  the  cylinder  is  given  as 


_  ‘  oo  f 

G(p,4')=-^  Y,  »«'(*«(>>)  p«(»^o/><)  +  ^(n)4*'(W<)]e--’”'*  (2.26) 


n=:— OO 


The  field  inside  the  coating  material  is  not  of  interest  in  this  research.  However,  in 
order  to  find  the  unknown  coefficient  .<4,(n),  it  is  necessary  to  find  the  field  inside 
the  coating  material  also.  Inside  the  coating  material  (i.c.  a.  <  p  <  b)^  G{p,<f>) 
should  satisfy 


1^  j9  J_  _£L  1.2 

P  dp  ^  dp  ^  p^  d4>^  ^ 


G(p,<i>)  =  0 


(2.27) 


It  is  assumed  that  the  source  is  located  outside  the  cylinder  (i.e.  p'  >  b).  Thus,  the 
total  field  G(p,<j>)  inside  the  coating  material  (i.e.  a  <  p  <  b)  can  be  expressed 
as 


(?(p,^)  = 


oo 

E 

n=:-O0 


Hl?\kop') 


B,in)Jn{kip)  +  D,{n)Hi^\kip)^  (2.28) 


The  unknown  coefficients  A,(n),  H,(n)  and  Dg{n)  can  be  determined  from  the 
boundary  conditions  given  as  follows; 


Ez{p,<(>)  =  0 


at  p  =  a  (2.29) 


Ez{p.<f>) 


is  continuous  at  p  =  b 


(2.30) 


is  continuous  at  p  s;  b 


(2.31) 


Thus,  the  scattering  coefficient  j4«(n)  is  given  as 


^,(n)  =  - 


Jn(kob)  -  jC,{n)Jn{kob) 
H^nhkob)-jC,{n)}ii^\koh) 


(2.32) 


•where  C,(n)  for  the  coated  cylinder  is  defined  as 


C(n)  =  -  -^  -  Nnik.aMkib) 

Jn{kia)Nn{kib)  -  Nn{kia)Mkib) 


(2.33) 


Zy  and  ki  denote 'the  characteristic  impedance  and  the  wave  number  of  the  coating 
material  respectively  (i.e.  Zy  =  =  ). 

It  is  interesting  to  note  that 


lim  C,(n)  =  CO 
a~-b 


(2.34) 


Equation  (2.34)  can  be  easily  proved  using  the  Wronskian  of  the  Bessel  functions. 
Thus,  when  the  tliickness  of  the  coating  is  very  small  in  terms  of  the  wavelength,  the 
behavior  of  the  field  is  essentially  the  same  as  for  the  TMt  case  of  the  conducting 
cylinder. 

By  replacing  Jn(a:)  and  Nn{z)  with  /fi^^(x)  and  Hn\x),  the  Equation  (2.26) 
can  be  rewritten  as 

+  (2.35) 


where 


F,(n)  =  - 


Hj^Ukob)  -  jC,{n)Hi%ob) 
H^'ikob)  -  jC,{n)Hi^\kob) 


(2.36) 
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Ct{n)  in  Equation  (2.33)  can  also  be  rewritten  ais 


]  rj  '  /ll  <21  (2)  /I)  (2.3  (  ) 

4  ^(A,6)^r(A:,a)  -  Hi:\kib)Hl^\kia) 

The  total  field  at  P{p  >  6,  (^)  for  an  incident  plane  wave  can  also  be  obtaiiied  from 
Equation  (2.35). 

Gip,4>)  =  ^  E  \Hi^\kop)-\-F,{n)H'^\kop)]  (2.38) 

^  n=— oo  ^  ■' 

2.3.2  TEx  case 

The  boundary  conditions  of  the  coated  cylinder  for  the  TEx  c^se  are  given  as 
follows; 

E^{p^4>)  =  Q  at  p  =  a  (2.39) 

ff2{pt4>)  *■»  continuous  at  p  —  b  (2.40) 

F^{p,4>)  *5  continuous  at  p  =  h  (2.41) 

By  using  the  same  procedure  as  for  the  TMx  case  with  the  above  boundary  con¬ 

ditions,  total  field  G{p,<f>)  is  given  as 

^  n=:— OO  ^  ^ 

(2.42) 


where 


=  - 


JUkob)-jCh{n)Jn{kob) 
Hi^hkob)  -  jC\{n)Hl?\kob) 


where  Cft(n)  for  the  coated  cylinder  is  defined  as 


=  '^nikia)Kikib)  -  KikMikib) 

^  Zo‘  JUkia)Nnikib)-  Ni,{kia)Jn{kib) 


It  can  be  easily  proved  that 


lim  C;,(n)  =  0. 
a— *5 


(2.43) 


(2.44) 


(2.45) 


Thus,  for  very  thin  coatings,  the  behavior  of  the  field  is  essentially  the  same  as  for 
the  TEg  case  of  the  conducting  cylinder. 

In  terms  of  the  Haiikel  functions  only,  the  Equation  (2.42),  (2.43),  and  (2.44)  are 
rewritten  as 


00  r 

=  £  rfW>)  +  (2.46) 

®  n=-co  ■* 


where 


Hi^^\kob)~jCh{n)Hi^\kob) 
Hl?^\kob)  -  jCh{n)Hi?\kab) 


(2.47) 


and 


(t,a) 

The  total  field  at  P(p  >  6,  ^)  for  an  incident  plane  wave  is  given  as 


(2.48) 


«((>.<«  =  ^  E  W''(i,f>)  +  F*(n)Hi^>{lb„(,)].-j“W-f)  (2,49) 

^  n=—oo  ■* 
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2.4  Examples  and  the  discussion 


Some  examples  of  the  eigenfunction  solution^  for  plane  wave  incidence  are 
given  in  Figure  2.1  ~  2.21.  The  geometry  of  the  problem  is  shown  in  Figure  2.1. 
Figure  2.2  ~  2.9  illustrates  how  the  normalized  field  pattern  of  the  impedance 
cylinder  changes  as  the  surface  impedance  increaises  from  zero  (i.e.  TEz  case)  to 
infinity  (i.e.  TMg  case).  In  these  examples,  surface  impedances  are  taken  to  be 
purely  imaginary  (i.e.  pure  capacitive  for  TAf*  smd  pure  inductive  for  TEz).  The 
cases  of  complex  impedance  will  be  given  in  later  Chapters  in  the  comparision  with 
the  ray  solution.  Figure  2.10  ~  2.15  and  Figure  2.16  ~  2.21  ilustrate  the  change 
of  the  scattering  pattern  as  the  thickness  of  coating  increases  for  TMz  and  TEz 
respectively.  In  these  examples,  €r=4  and  pr=l-  The  cases  of  complex  tr  and  pr 
are  given  in  later  Chapters. 

It  ran  be  observed  that  manifest  effects  due  to  the  changing  values  of  the  impedance 
or  the  thickness  are  shown  for  the  field  pattern  in  the  shadow  region.  From  these 
examples,  we  can  observe  that  the  scattered  field  in  the  shadow  region  does  not 
simply  increases  or  decreases  as  the  surface  impedance  or  the  thickness  of  the 
coating  increases.  However,  it  is  not  possible  to  give  any  physical  intepretation 
for  the  change  of  the  scattering  pattern  from  the  eigenfunction  solutions  before  we 
introduce  the  ray  concept  which  will  be  given  in  later  Chapters.  Thus,  detailed 
discussion  of  the  scattering  mechanism  is  deferred  to  later  Chapters. 


*The  FORTRAN  subroutine  for  the  Bessel  functions  provided  by  Prof.  J.H.  Richmond  have  been 
essential  for  obtaining  these  results. 
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Figure  2.3:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 


Figure  2.4:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 
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Figure  2.6:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 


Figure  2.7:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 
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Figure  2.8:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 
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Figure  2.9:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 


Figure  2.10;  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TMt  case 


Figure  2,11:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TMz  case 
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Figure  2.12:  Normalized  bistatic  scattering  pattern  of  coaled  cylinder  :  TM- 


Figure  2.13:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TMz 
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Figure  2.14:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TM-  case 
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Figure  2.15:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TMz  case 
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Figure  2.16:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TE:  case 
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Figure  2.17:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEt  case 
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Figure  2.18;  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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Figure  2.19:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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*o*  =  20  «<  =  0.5A<, 


Figure  2.20:  Normalized  bistatic  scattering  pattern  of  coaled  cylinder  :  TEz  case 
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CHAPTER  III 


GTD  SOLUTION 


Ray*optical  solutions  for  the  Held  scattered  by  the  impedance  and  the  coated 
cylinder  are  obtained  in  this  Chapter.  For  the  perfectly  conducting  cylinder,  Keller 
[1]  developed  an  asymptotic  solution  to  this  problem  and  interpreted  it  in  the  ray- 
optical  terms  within  the  framework  of  his  geometrical  theory  of  diffraction  (GTD). 
Even  though  the  boundary  conditions  to  be  considered  are  different,  the  same  pro¬ 
cedure  as  that  for  the  derivation  of  GTD  solution  for  the  conducting  cylinders  can 
be  taken  to  deduce  the  GTD  solutions  for  the  impedance  and  the  costed  cylin¬ 
der.  First,  the  Held  representation  in  the  angular  eigenfunction  series  as  given  in 
Chapter  II  is  transformed  into  the  integral  form.  Then,  the  integral  is  evaluated 
asymptotically  vis  the  stationary  phase  method  or  the  residue  theorem. 

It  is  noted  that  the  behavior  of  high-frequency  field  is  characterized  by  various 
geometric  optical  domains  as  illustrated  in  Figure  3.1.  Separate  analysis  is  given 
for  each  geometric  optical  domain  as  the  field  representation  has  different  format 
according  to  the  region  where  the  field  point  is  located.  Lit  region(I)  and  shadow 
region(III)  are  separated  by  the  shadow  boundaries  (SB)  according  to  the  illumi¬ 
nation  from  the  source.  Small  angular  ranges  adjacent  to  the  shadow  boundaries 
which  are  indicated  as  shaded  region(II)  are  the  prenumbral  or  transition  region. 
The  solution  obtained  in  this  Chapter  is  not  valid  in  this  region.  The  analysis  for 
the  field  in  transition  region  will  be  given  in  Chapter  V.  In  this  analysis,  we  also 
exclude  the  surface  or  caustic  boundary  layer  region(IV)  which  is  in  the  immediate 
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vicinity  of  the  cylinder  surface. 

As  mentioned,  Keller’s  pure  ray-optical  solution  fails  within  the  transition  region 
adjacent  to  the  shadow  boundaries.  The  failure  of  Keller’s  GTD  solution  has  been 
overcome  by  the  development  of  a  uniform  version  of  GTD  by  Pathak  (3).  Hence, 
it  may  seem  to  be  obsolete  to  develop  the  pure  ray-optical  solution  for  our  problem 
which  will  inevitably  fail  in  the  transition  region.  Nevertheless,  it  appears  that  the 
ordinary  GTD  solution  is  a  necessary  step  to  be  taken  for  the  development  of  the 
uniform  GTD  which  will  be  discussed  in  Chapter  V.  This  is  because  the  ray  format 
to  be  used  in  the  uniform  GTD  solution  comes  from  the  ordinary  GTD  solution. 
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REGION  (S) 


TRANSITION 


3.1  Integral  Representation  of  the  Field 

As  discussed  in  Appendix  A,  the  field  representation  in  the  eigenfunction 
series  form  given  in  Chapter  II  can  be  transformed  into  an  integral  form.  The 
eigenfuntion  series  solution  of  interest  has  the  form  given  by 

G{p,<t>)  =  u=  .B,^h(n)exp(-jTi^)  (3.1) 

n=— oo 

where  Bf  fi(n)  depends  on  the  polarization  of  the  incident  field  and  the  subscripts 
s  and  h  corrrespond  to  the  TM^  and  TBz  cases  respectively. 

From  Equations  (2.13),  (2.19),  (2.35)  and  (2.46),  B,^f,(n)  is  given  as 


B,,k(n)  =  Hi^\koP>) 


(3-2) 


From  Equation  (A.3)  in  Appendix  A,  the  eigenfunction  series  can  be  transformed 
into  the  integral  representation. 

H  exp(-jn^)  =  B,  h{u)  ---  ^  dv  (3.3) 

n=-oo  2  JCt  '  Sini/TT 

where  C*  is  the  contour  which  encloses  the  lower  half  of  complex  r^-plane  with  a 
large  semicircle.  From  Equation  (3.2),  Bg  f^{u)  is  given  by 

B.^f)  =  Y  [»i‘’(W<)  +  F,,kW)  (3.4) 


and 


(3.6) 
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For  the  impedance  cylinder, 


C,{u)  =  C,  =  g 


far  TMt 


Chiu)  =  Ch^^ 


far  TEt 


For  the  coated  cylinder,  C«(t/)  and  Ch(u)  &re  ipven  as 


Cs(u)  = 


far  TMj 


Ch{u)  3  ^  /i\  "  (oy  {2\  (iV  I  Ex 

HV\kib)Hir^  (*lfl)  -  Hi>\kib)ffy^  {kia) 


It  is  noted  that  ^i^j(x)  =  and  /rS(x)  =  e">*"^i^^(x) ,  there* 

fore,  it  can  be  easily  seen  that  is  an  even  function  of  u  (i.e.  i/)  = 

From  the  Equation  (A.5)  in  Appendix  A, 


/QO-J< 

.  S$,hi'') 

-oo~j( 


cos  u[<f>  —  x) 
sini/x 


(3.10) 


co$u{^  —  x) 
sint/x 
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It  is  noted  that 


(3.11) 


cos  i/(^  -  tr)  .  g->i/(2ir+^)  ^  ^-ju(2w~^) 

sinj/TT  ~  l_e-j2t/)r 


(3.12) 


Thus,  the  integral  representation  of  the  field  can  be  decomposed  into  two  terms  as 


O  I  cw 

U  =  tt  +  u 


(3.13) 


where 


^-M2k+<P)  g-ji/{2jr-^) 

X  - rs -  dly 

1  _  e-j2i/ir 


(3.15) 


It  has  been  pointed  out  [13]  that  u®  in  Equation  (3.14)  contains  all  the  geometric- 
optical,  transition  and  the  dominant  dififraction  effects.  In  the  following  sections, 
the  integral  in  Equation  (3.14)  will  be  evaluated  asymptotically  via  the  stationary 
phase  method  to  yield  the  geometric-optical  fields  in  the  lit  region.  A  residue 
series  solution  for  the  integral  gives  the  creeping- wave  representation  in  the  shadow 
region.  The  lit  and  the  shadow  region  separated  by  the  shadow  boundary  (S.B.) 
are  indicated  in  Figure  3.1.  The  term  u®®  represents  multiply  encircling  creeping 
waves,  which  will  be  discussed  later  in  this  Chapter. 
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3.2  Field  in  the  deep  lit  region  >  G.O.  field 

Lit  region  is  in  the  angular  range  of  0  <  \^\  <  ,  where  the  field  point 

is  in  the  line  of  sight  from  the  source.  In  this  region,  the  asymptotic  evaluation 
of  given  in  Equation  (3.14)  leads  to  geomeric-optical  fields;  incident  field  «* 
and  reflected  field  vT .  It  should  be  noted  that  this  solution  is  valid  in  the  deep  lit 
region  excluding  the  transion  region  in  the  vicinity  of  the  shadow  boundary.  For 
simplicity  of  calculation,  we  assume  that  p  <  p' .  For  the  case  of  p  >  p\  the  same 
procedure  can  be  taken  to  obtain  the  same  result. 

Equation  (3.14)  is  decomposed  into  two  terms  as 

u°  =  uf-|-uf  (3,16) 

where 

uf  =  H^^\kop')Hi^\kop)e-j''^du 

O  J-CC-Jt 

and 

o  J-oo-jt 

3.2.1  Evaluation  of  u| 

Equation  (3.17)  can  be  asymptotically  evaluated  via  the  stationary  phase 
method  after  H^\kop')  and  Hi}\kop)  are  replaced  by  their  asymptotic  ap¬ 
proximations.  The  problem  with  this  method  is  that  the  Hankel  functions  have 
difierent  asymptotic  approximations  in  the  different  regions  along  the  path  of  the 
integral  where  u  changes  from  -00  to  -f  00.  Therefore,  strictly  speaking,  it  is  neces¬ 
sary  to  devide  the  integral  path  into  several  sections  according  to  the  valid  regions 
of  the  approximation,  which  makes  the  evaluation  cumbersome.  To  circumvent 


(3.17) 


(3.18) 
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this  difficulty,  an  alternative  method  will  be  employed  to  obtain  the  same  result. 
First,  we  represent  Hj^\kop')  and  Hl^\kop)  by  their  integraJ  form. 


H^\kop')  =  ~  /  exp(— jA:p*sina  -f  jva)  da 
rr  yCj 


(3.19) 


H^\kop)  =  —  f  t:xp{- jkp sin  13  +  jV/?)  d0  (3.20) 

K  JCj 

The  contour  C\  and  C2  are  depicted  in  Figure  3.2.  After  replacing  H^\kop*)  and 
H\}\kop)  in  Equation  (3.17)  by  their  integral  forms  and  rearranging  the  integral, 


we  obtain 


“1  =  ^  -  jkop'  sin  a) 


^  -  ii>)]du\  d(3  da 


(3.21) 


Assuming  that  a  and  /3  are  real  variables,  it  can  be  seen  from  the  well-known 
relationship  of  the  inverse  Laplace  transformation  that 

/00-jf  /■«+><» 

.  exp[jV(Q  +  /3-  <^)]  di/=  -j  /  exp[s(o(  +  0  -  <f>)]  da 
-00-jt  Je—joo 


=  2Tt8{a  -f- 


Incorporating  Equation  (3.22)  into  Equation  (3.21),  we  obtain 


(3.22) 


Uj  =  ^  ^  exp[-jfco/>^sina  -  jkopsm{<i>  -  a)]  da 


(3.23) 


In  obtaining  Equation  (3.23),  we  assumed  that  <(>  —  a  lies  on  the  contour  Cj. 
Thus, 

-ir<^-o<0  (3.24) 


Asymptotic  evaluation  of  the  integral  in  Equation  (3.23)  can  be  easily  obtained 
via  the  stationary  phase  method.  Equation  (3.23)  can  be  rewritten  as 


l^Ici  W'^(a)]  da 


(3.25) 


where 


#(Qr)  =  —  sin  o  —  ^  sin(^  —  a) 


(3.26) 


At  the  stationary  point,  ^^(a«)  =  0,  therefore,  the  stationary  point  a,  is  deter¬ 
mined  by  the  relationship 


$'(as)  = 


a#(Q) 


da 


= -cosa,  + -f  cps{4>  ~  a,)  =  0.  (3.27) 

aszQg  P 


The  stationary  point  should  lie  on  the  contour  C2  in  order  for  u  j  to  have  nonzero 
value.  Therefore 


0  <  a,  <  TT  (3.28) 

From  Equations  (3.24)  and  (3.28) 

4>  <  a,  <  It  (3.29) 

It  was  assumed  p'  >  p  and  0  <  ^  <  tt.  Thus,  it  can  be  shown  that  in  order 
for  OTg  to  satisfy  the  given  condition  in  Equation  (3.29) 
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meaning  of  the  stationary  point 


*  <  CO.  ■ 


(3.30) 


Thus,  Uj  contributes  only  when  the  above  condition  is  satisfied. 
Final  result  of  the  asymptotic  evaluation  of  the  integral  yields 


(3.31) 


The  physical  meaning  of  the  stationary  point  is  illustrated  in  Figure  3.3.  From  the 
Figure  3.3,  we  observe  that  the  distance  from  the  source  to  field  point  is  ^ven  as 


I  =  \p'  -  ^  p  sma,  -  p sin(Qr,  -  <i>) 


(3.32) 


Therefore,  u®  can  be  expressed  as 


e-jj 

2v/5^  y/iej 


for  0  <  (ff  <  cos' 


■'&) 


(3.33) 


From  the  above  result,  we  conclude  that  is  the  direct  incident  field  from  the 
source  to  the  field  point  for  0  <  ^  <  cos~^{p/p^).  In  obtaining  above  result, 
the  integral  representations  for  the  Hankel  functions  are  used  instead  of  the  Debye 
approximations.  Even  though  the  same  result  can  be  obtained  by  using  the  Debye 
approximations,  this  method  is  mathematically  more  rigorous. 

Note  that  u°  is  independent  of  the  boundary  conditions  of  the  cylinder.  Thus  the 
result  obtained  here  can  be  applied  to  both  the  impedance  and  the  coated  cylinder. 
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S.2.2  Evaluation  of 

In  this  section,  112  Equation  (3.18)  is  evaluated  asymptotically  via  the  sta¬ 
tionary  phase  method.  The  alternative  method  used  for  the  evaluation  of  u|  can 
not  be  easily  applied  here  because  in  the  integrand  is  a  function  of  u. 

Therefore,  we  replace  the  Hankel  functions  in  the  integrand  by  their  asymptotic 
approximations.  It  is  noted  that  the  Hankel  functions  in  the  integrand  have  dif¬ 
ferent  asymptotic  approximations  accciding  to  the  changing  vsdues  of  v  along  the 
contour  path.  In  order  to  take  a  mathematically  rigorous  procedure  it  is  necessary 
to  divide  the  integral  path  into  many  sections  because  the  integrand  of  Equation 
(3.18)  has  Hankel  functions  with  3  different  arguments  (  A:©/?,  A^oP^  kc,b  )  for  the 
impedance  cylinder  and  the  5  different  arguments  (  koP,  kop\  kob,  kia,  kib) 
for  the  coated  cylinder.  To  simplify  the  procedure  we  first  divide  the  integral  into 
two  terms  as 

“2  =  +  ^2  (3.34) 

where 

o  Jkob-j(  ' 

and 

O  J-O0-}€ 

Evaluation  of  Ii 

First,  we  assume  that  the  asymptotic  evaluation  of  the  integral  depends  only 
on  the  stationary  phase  point  and  therefore  the  contribution  from  the  end  point 
at  1/  =  kob  can  be  ignored.  This  assumption  is  valid  when  the  stationary  phase 


(3.35) 


(3.36) 
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point  v$  is  not  close  to  koh.  When  the  field  point  is  in  the  deep  lit  region,  it  can  be 
shown  that  the  stationary  phase  point  is  located  far  enough  from  koh.  However, 
as  the  field  point  approaches  the  shadow  boundary,  this  assumption  is  no  longer 
valid.  The  case  when  the  approaches  koh  will  be  discussed  m  Chapter  V  in  the 
treatment  of  the  field  in  transition  region. 

Assuming  that  v  is  not  close  to  kjb  in  the  given  integral  path  of  Equation  (3.35) 
where  koh  <  Re[v\  <  oo,  it  can  be  observed  that 

ffl^hkoh)  % -ffi^\kob)  (3.37) 


H^u\koP)H]?\kop)e~^^*du  (3.39) 

8  J  kob~jf 

The  integrand  in  above  Equation  (3.39)  includes  two  Hankel  functions,  Hl^\kop') 
and  Hi/  '{kop).  As  remarked  before,  they  have  different  asymptotic  approximations 
along  the  path  of  the  integral.  Here,  we  make  another  assumption  that  both  p’ 
and  p  are  large  enough  that  the  stationary  phase  point  u,  is  always  located  such 
that; 

<<  1  and  7^  <<  1  (3.40) 

kop 

This  is  a  reasonable  assumption  because  the  field  in  the  surface  boundary  layer 
region(IV)  of  Figure  3.1  is  not  of  interest  in  this  analysis.  From  the  Debye  approx¬ 
imations  of  the  Hankel  functions  for  this  range  of  u 
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(3.41) 


Jll%oP')d^\koP) 


2j  exp[-jA:ojO%in  a  —  jkopt,m0  +  juja  -f  0)] 
yjp'p  siua  #in/? 


rk, 


where  a  =  cos  ^{t^/kop')  and  ^  ~  cotr^(u/kop). 
Therefore,  /j  can  be  approximated  as 


h 


1  f  exp[jtop'»(v)] 
Airko  JCy  sino  sin/? 


(3.42) 


where  ^{v)  is  defined  as 


$(i/)  =  —  sino  — 


4sin/3  +  <i>) 

P  P 


(3.43) 


and  Cy  denotes  the  integral  path  of  v  subject  to  Equation  (3.40)  and  the  condition 
that  \vlkoh\  »  1.  At  the  stationary  phase  point,  =  0.  Therefore,  it  can 

be  easily  shown  that 

oi$+0$  =  4>  (3.44) 


where  o,  =  cos~^{i/g/kop')  and  fig  —  cos“*(i/j/AoP)-  Physical  meaning  of  the 
Q(  and  fig  is  depicted  in  Figure  3.4.  It  is  important  to  note  that  both  <Xg  and  0g 
should  satisfy  the  condition; 


0  <  i2c[a,]  <  TT  and  0  <  Jtel^g]  <  ir  (3.45) 

It  can  be  verified  that  in  order  for  a,  and  /3g  to  He  on  the  integral  path  Cy  subject 
to  the  above  condition,  the  following  condition  on  (f>  should  be  satisfied 

cos-1  ^  ^3  46) 
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p 


Figure  3.4:  Physical  meaning  of  the  stationary  point 
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where 


4>,h  -  cos' 


■■(?) 


+  cos' 


(3.47) 


For  the  region  where  <f>  satisfies  the  above  condition,  the  asymptotic  evaluation  of 
the  integral  in  Equation  (3.42)  can  be  obtained  via  the  stationary  phase  method. 
The  result  is  given  by 


cxpl-j{kop*  sin  a$  -f-  kpp  sin  $,)] 
yjkop'  »\na$  +  koptinfit 


(3.48) 


From  the  Figure  3.4,  it  can  be  seen  that 


^  =  p'siuQ,  +  psinfit 


(3.49) 


Therefore, 

e~^  4  ^  (  p\ 

'TU  b  j  ^ 

Thus,  it  is  shown  that  Ii  gives  the  same  result  as  the  evaluation  of  Uj  for 
0  ^  ^  <  cos“^(p/p') . 

3.2.3  Incident  field 

We  conclude  from  the  results  we  obtained  so  far  that  the  incident  field  u* 
consists  of  two  terms.  For  0  <  ^  <  cos~^(p/p*),  u*  comes  from  the  asymptotic 
evaluation  of  ul  given  in  Equation  (3.33)  and  for  co8“*(p/p')  <  ^  <  ^,i>,  tt* 
comes  from  the  evaluation  of  /j  in  Equation  (3.48).  Consequently,  we  obtain  the 
incident  field  «,•  for  the  lit  region  (0  <  ^  <  ^,5)  as 


g-jf  g-jkd 

2V^  \/k^ 


for  0  <  <l>  < 


(3.51) 
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This  result  is  valid  for  both  the  impedance  and  the  coated  cylinder.  In  fact,  the 
direct  incident  from  the  source  is  independant  of  the  boundary  conditions  of  the 
scatterer.  It  should  be  emphasized  again  that  these  results  are  not  valid  when 
the  field  point  is  close  to  either  the  shadow  boundary  (S.B.)  or  the  surface  of 
the  cylinder  where  the  Debye  approximations  of  the  Hankel  functions  cannot  be 
applied. 

3.2.4  Evaluation  of  I2 

A  similar  procedure  as  used  for  the  evaluation  of  I\  can  be  applied  to  obtain 
the  asymptotic  solution  for  I2.  First  the  Hankel  functions  in  the  integrand  are 
replaced  by  their  appropriate  asymptotic  approximations  and  then  the  stationary 
phase  method  is  used  to  obtain  the  final  result. 


/2  =  -/  (3.52) 

o  J-ao-jf 


Equation  (3.5)  can  be  rearranged  as 


Fs,h{^)  = 


(  d'Uw/d^hkob)-jc,,h{^) 

\  Hl^^\kob)/d^\kob)  -  jC,^h{^) ) 


(3.53) 


We  define  Sg  h{u)  as 


—  (2)1  (21 

nr  {kob)/Hl^\kob)  -  jCg^hH 


(3.54) 


From  the  Debye  approximations  of  the  Hankel  functions,  it  can  be  shown  that  for 
—00  <  1/  <  kob, 


dJ\kob) 


irkob  sin  7 


exp  |j  ^feobsin  j  -  ji/  -  (3.55) 
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rU) 


H^ikob)  .  .  . . 

— TXv - ss  -;  exp{2j(M*m7  -  71/)] 


M  .  . 

-711 - «s  J«n7 

Hy\kob) 

and 

— T2X - %  -jsin7 

jyr^ibo6) 

where  7  =  cos~ {vjkob). 

Thus,  Fg^fi{u)  and  5,^h(i/)  can  be  sinsplified  as 


^t,hi^)  ^  -3  exp[2i(Msm7  -  71/)] .  5,,ft(i/) 


where 

e  /  \  sin 7  -  C,,h(i/) 

*i«7  +  C',,fc(i/) 

Substituting  Equation  (3.60)  into  Equation  (3.52),  we  obtain 

/2  a  ^  4*>(»:.p')ffi**(*o,>)«p[2j(M»in7  -  7-')  - 

O  «/— 00— J€ 


By  using  the  Debye  approximations  for  H^\kop*)  and  Hi^\kop)^  I2 
pressed  as 


(3.57) 

(3.58) 

(3.59) 

(3.60) 

(3.61) 

-  jV^]  du 

(3.62) 
can  be  ex- 
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-j  fkob-jt  S,^hii')  ^ 

^2  =  TZT  /  .  r. . : - -  ;  cxp\j  hop  $(»/))  du 

^nko  J-oo~}t  ^ p  sin  a  Bin  0 


(3.63) 


where 


p  26  u 

$(i/)  =  —  sin  Q - r  sin/3  +  — T  sin  7  +  - — (a  +  0  —  2'f  —  4^) 

P  P  KoP 


(3.64) 


The  derivatives  of  ^(t')  with  respect  to  */,  are  given  as 


(3.65) 


=  J_  f  _J _ 1 _ L_) 

di/2  kop'  \  fco6  sin  7  kop^  sin  a  kop  Bin0  / 


(3.66) 


In  order  to  obtain  the  derivatives^  we  used  the  relationship 


u  =  kop'  cos  a  ^  kop  cos  0  =  kob  cos  7 


(3.67) 


therefore 


do  _  —1 

du  koP'  sin  o 


ii.  - 

du  koP  sin  0 


dj  _  —1 

du  kob  sin -f 


(3.68) 


The  stationary  phase  point  u,  is  determined  by  the  relationship  given  by 


aa  +  0s-  27,  =  <f> 


(3.69) 


where  o,  =  cos  ^{i^g/kop'),  0a  =  cos  ^{u,/kop)  and  7,  =  cos"*(i/*/A:o6), 
At  the  statioirary  point  Ua 
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(3.70) 


2b 

^{y$)  —  ^*in7*  —  sin  a. 


p' 


and 


The  physical  meaning  of  the  stationary  phase  point  is 
should  be  noted  that 


kob  sin  ^0  kop*  sin  a0 


Bin  a) 


kop  sin 


(3.71) 


depicted  in  Figure  3.5.  It 


0  <  Jic(a,]  <?r,  0  <  Re[$0]  <  ir,  0  <  i2e|7,]  <  ir  (3.72) 


From  the  Figure  3.5,  it  can  be  seen  that 


p'sinctf  >  &sin7« 

psin^s  >  5  sin  7,  (3.73) 


Therefore 


Viv0) 


p'  sin  a,  —  fesin7j  ^  pan  fig  —  58in7» 
kop'b  sin  ag  sin  7*  k^pb  sin  fig  sixt-^g 


(3.74) 


The  asymptotic  evaluation  of  the  integral  using  the  stationary  phase  method  leads 
to  the  result 


h 


zi 

4 


Sg,hM 


exp(i(2jbo6sin7«  —  kgpanfig  —  X;op'sina«)] 


X 


/ 


68in7j 

Ikop'p  sin  <xg  sin  fig  —  kopb  sin  fig  sin  7^  —  kop'b  sin  ag  sin  jg 


After  the  rearrangement  of  the  above  equation,  we  obtain 


(3.75) 
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^2  ~  X  V »  - /.  ,  .  ...  - 

’  »  t  /JboP  sin  a*  —  *©6  sm  7, 


frsm7«(p'8ina«  —  68in7«) 


2(p'  sin  a,  —  6  sin  7j  )(p  sin  —  5  sin  7j  )  +  fc  sin  sin  ct$-¥  p  sin  —  2b  sin  7*  ) 


X  €Xp[—j{kop^n0t  —  A;ofrsin7«)] 

FVom  the  Figure  3.5,  we  observe  that 


(3.76) 


ll  =  p'  sin  at  —  b  sin  7* 


(3.77) 


I2  =  psiTifi,  —  56in7» 


(3.78) 


Thus,  1%  can  be  expressed  as 


r  _  e-iMi  ^  ^  I  12b  sin  ft 

2v^  VW  ’  V  2/,/2  +  (ii  +  /2)6«n7,  * 


(3.79) 
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3.2.5  Reflected  Field 


The  result  of  the  asymptotic  evaluation  of  I2  as  given  in  Equation  (3.79)  rep¬ 
resents  the  field  reflected  on  the  surface  of  the  cylinder.  From  the  result  obtained 
in  the  previous  section,  it  is  clear  that 


2v/2ff 


=  u\Q) 


(3.80) 


where  «*(Q)  represents  the  incident  field  at  the  specular  reflection  point  Q  as 
shown  in  Figure  3.6.  From  the  comparision  of  this  result  with  the  ray  format 
of  the  perfectly  conducting  cylinder  it  is  apparent  that  i5,,fc(i'j)  represents  the 
reflection  coefficient  Rg^h.' 

Thus  Equation  (3.79)  can  be  rewritten  as 


/  y  (3.81) 

2v^  vC?r  ’  V  2/1/2  +  (/i  +  h)bcos  Bi  ' 


We  define  Sp  as 


Sp  = 


I  lib  cos  Bi 

2/1/2  +  (/l  +  /2)^cos  Bi 


(3.82) 


It  can  be  easily  identified  that  Sp  is  the  usual  spreading  factor  associated  with  the 
reflection  field  from  a  curved  surface,  which  is  given  as 


(3.83) 


where  the  reflected  ray  caustic  distance  p’’  can  be  calculated  from  the  relationship, 


J-  _  i  2 

p’’  /i  6  cos  Bi 


(3.84) 
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The  spreading  factor  Sp  accounts  for  the  spreading  of  electromagnetic  energy  after 
the  reflection  occurs. 

Thus,  the  ray  format  for  the  reflected  field  can  be  expressed  as 


u 


=  «■(<?)  •  li..h  ■  Sp 


(3.85) 


It  is  interesting  to  observe  that  the  ray  format  for  the  reflected  field  associated 
with  the  impedance  or  the  coated  cylinder  remains  the  same  as  that  associated 
with  the  perfectly  conducting  cylinder. 

Reflection  coefficient 

Since  the  reflected  ray  should  satisfy  Snell's  law  of  reflection,  the  incident 
angle  $i  must  be  equal  to  the  reflection  angle  6r  nt  the  point  Q  in  the  Figure  3.6. 
From  the  physical  meaning  of  the  stationary  point  7,,  it  is  noted  that  7,  :s  ^  ~ 
(i.e,  sin 75  =  cos^,).  Thus,  the  reflection  coefficient  for  the  impedance  cylinder  is 


cosej- C,^h 

cos^j  +  C,,;, 


(3.86) 


Substituting  the  Equations  (3.6)  and  (3.7)  into  Equation  (3.86),  the  reflection 
coefficients  for  the  impedance  cylinder  are  given  as 


Zq  -  Zg  cos  Sj 
Zo  +  Zf  cos  6i 


far  TMt  (3.87) 


and 


rt  -  Z,/cos6i 

Zo-\-Z,/cos9i 


far  TEz  (3.88) 


Note  that  the  reflection  coefficients  of  the  impedance  cylinder  are  the  same  as  the 
well-known  reflection  coefficients  for  a  plane  wave  incidence  on  a  infinite  ground 
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0 


(a)  Impedance  Cylinder 

P 


(b)  Coated  Cylinder 


Figure  3.6:  Reflection  field 


plane  with  surface  impedance  Z«. 


As  ^ven  in  the  Appendix  B,  the  reflection  coefficients  for  the  coated  cylinder 
can  be  obtained  by  the  asymptotic  approximation  for  Sg^h{ut)  using  the  Debye 
approximations  of  the  Hankel  functions.  The  results  are  summarized  as  follows 


Zq  —  j Zi  cos Bj / cos tan{^’(a» Kki)} 
Zo  +  jZi  cos  0i/ cos  $i  tan{4'(a»6,fcl)} 


for  TM  (3.89) 


and 


_  _  Zq  ~  jZi  cos  0if  cos  Oj 
^  +  jZi  cos  0i/  cos 6i  tan{'i'(a, 6,  kj )} 


for  TE  (3.90) 


where  Z\  represents  the  characteristic  impedance  of  the  coating  material  (i.e. 
Zi  —  y/S  j  and  and  St  are  the  incident  and  the  transmission  angle  of  the 
wave  at  the  air/coating  interface.  '9{a,b,ki)  is  defined  as 


9{a^b,k\)  =  fci(6sin/?fc~  asin/fa)  -  kob[0i,  -)0o)8in^i  (3.91) 


where 


and 


(3.92) 


(3.93) 


Comparing  the  reflection  coefficients  for  the  coated  cylinder  with  those  for  the 
impedance  cylinder  it  can  be  observed  that  the  equivalent  impedeuice  of  the  coated 
cylinder  associated  with  the  reflected  field  is 
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=  iZi  ■ 

COS^t 


for  TM^ 


(3.94) 


cos  ^<tan{($(a, 6,^1 )}  for  TEz 


(3.95) 


This  result  shows  that  the  equivalent  impedance  is  not  only  dependent  on  the 
polarization  of  the  incident  field  but  also  on  the  incident  angle  at  the  reflection 
point.  Therefore,  it  is  apparent  that  the  constant  surface  impedance  cylinder  can 
not  be  appropriately  used  for  the  equivalent  model  of  the  coated  cylinder.  However, 
it  will  be  shown  later  in  this  Chapter  that  the  equivalent  impedance  in  association 
with  the  diffracted  field  is  constant  along  the  cylinder  surface.  For  the  thin  coating, 
it  can  be  seen  that  0a  ^  0b  therefore  $(0, 6,  hj)  can  be  approximated  as 


'9{a,b,ki)  X  ki(b-  a)  sin  0a  =  kid  cos  0t 


(3.96) 


Thus,  for  a  thin  coating,  Rg  and  R^  can  be  approximated  as 


_  Zq  —  jZi  cos  Oj/ cos  $i  tan{kidcos$i) 
Zo  +  j Zi  cos  /  cos  0i  tan( fci  d  cos  Of ) 


for  TM  (3.97) 


Zq  —  jZi  cos  Qj/ cos  $i  tan(fe| d cos  9i ) 
Zo  +  jZi  cos  $t  /  cos  0i  tan(A;]  d  cos  6% ) 


for  TE 


(3.98) 


zi^W  =  jZi  ■ 

cos^< 


for  TMz 


(3.99) 


—  j^l  cosSi  t&n{kidcosBi) 


for  TEx 


(3.100) 


It  can  be  noticed  from  Equation  (3.97),  (3.98),  (3.99)  and  (3.100)  that  the  re¬ 
flection  coefficients  and  the  equivalent  surface  impedances  of  the  coated  cylinder 
coincide  with  those  of  the  grounded  planar  diclectric/ferrite  slab  with  plane- wave 
incidence  [7],  This  result  is  consistent  with  the  postulate  of  the  GTD  solution  that 
the  reflection  is  a  local  phenomenon  at  high  frequency. 

It  is  well-known  that  the  reflection  coefficients  for  the  conducting  surface  are 
fixed  to  ±  1  regardless  of  the  incident  angle  at  the  reflection  point  Qr.  However 
for  the  impedance  cylinder,  the  reflection  coefficients  are  dependent  on  the  incident 
angle  6^  and  the  surface  impedance  Z,.  For  the  coated  cylinder,  the  refection 
coefficients  are  dependent  on  the  incident  angle  as  well  as  on  <r,  /ir  and  thickness 
d  of  the  material  coating.  Figure  3.7  and  3.8  illustrate  the  change  of  of  the 
coated  cylinder  as  Cr?  Mr  and  d  changes  for  the  normal  incidence  (i.e,  6^  —  0). 
For  the  lossless  material  coating  where  both  (r  and  fir  have  real  values,  it  can  be 
easily  seen  that  the  phase  of  .R,  ^  changes  as  d  changes  but  =  1  regardless 

of  coating  thickness.  Physically,  the  phase  of  R,  accounts  for  the  lagging  of  the 
reflected  field.  If  there  exists  a  loss  factor  in  «r  or  fir  of  the  coating  material, 
can  be  smaller  than  1.  Consequently,  the  magnitude  and  the  phase  of  the  reflected 
field  can  be  controlled  by  proper  choice  of  the  relevant  parameters  (i.e.  e,.,  /i,.,  and 
coating  thickness  d). 
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I’'*'?/! 


ThicknesE  of  coaling 
tt-T  =  1- 

Cr  =  5.  -  >0  •“  ■ 

fr  —  5*  *  jO.l 

Cr  =  5.  -  iO.5  —  —  —  —  — 

Cr  =  5.  -  jl.  . . * 

m  rntimmmm  I 

Cr  =  5.  -  j5. 

Figure  3.7:  Reflection  Coefficient 
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Thickness  of  coating 


Figure  3.8:  Reflection  Coefiicient 


3.3  Field  in  the  deep  shadow  region 

In  this  section  a  ray  format  for  the  field  in  the  shadow  region  is  obtained  for 
both  the  impedance  and  the  coated  cylinder.  As  discussed  in  the  previous  section, 
the  behavior  of  the  geometric-optical  field  (i.e.  incident  field  and  reflected  field) 
for  the  impedance  and  the  coated  cylinder  are  almost  the  same  as  for  the  perfectly 
conducting  cylinder  except  for  the  difference  in  the  reflection  coefficients.  Manifest 
effects  of  the  surface  impedance  and  the  coating  of  the  dielectric/ferrite  materials 
can  be  observed  in  the  diffracted  field  in  the  shadow  region. 

The  diffracted  field  is  entirely  associated  with  the  creeping  wave  on  the  curved 
surface  of  the  cylinder.  The  creeping  wave  which  propagates  azimuthally  on  the 
cylinder  surface  carries  the  electromagnetic  energy  into  the  shadow  region.  A 
detailed  discussion  of  the  characteristics  of  creeping  wave  will  be  given  in  the  next 
Chapter. 

3.3.1  Residue  series  solution 

The  diffracted  field  can  be  obtained  by  a  residue  series  solution  for  the  integral 
representation  of  the  field  given  in  Equation  (3.14).  As  discussed  in  Appendix  A, 
the  residue  series  representation  of  the  field  is  the  result  of  a  Watson  tranformation 
of  the  field  given  in  eigenfunction  series  form.  As  illustrated  in  Figure  3.9,  the 
integrzd  path  of  the  Equation  (3.14)  is  deformed  such  that  the  contour  C*  encloses 
the  lower  half  of  complex  i/-plane  with  a  large  semicircle.  It  can  be  shown  that  the 
contribution  from  the  arc  of  the  large  semicircle  vanishes  as  its  radius  approaches 
infinity.  Thus,  one  obtains 

u'’  =  ^  /c.  e~^''^du  (3.101) 
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Ifn[v] 


Figure  3.9:  The  contour  of  the  integral 


The  integral  in  Equation  (3.101)  is  now  evaluated  using  the  residue  theorem.  Be¬ 
fore  we  use  the  residue  theorem,  the  poles  of  the  integrand  must  be  found.  The 
poles  of  the  integrand  come  only  from  the  poles  of  Fg  h(u).  From  the  Equation 
(3.5),  the  poles  of  Fg^h{u)  are  obtained  by  finding  the  roots  of 

=  0  (3.102) 


Equation  (3.102)  is  referred  to  as  transcendanial  equation.  A  detailed  discussion 
of  the  roots  of  the  transcendantal  equation  is  given  in  Chapter  IV.  The  residue 
series  solution  subject  to  Equation  (3.102)  is 


^  n=l 


d 


d]}\kcb)  -  ic,M'',)d]}(kob) 


(3.103) 

This  equation  can  be  further  simplified  by  using  the  Wronskian  of  the  Hankel 
functions.  From  the  relationship  of  the  Wronskian, 


W  =  dl\k<,b)«!^\Kb)  -  Bil>'(kob)nll>(kob)  =  ^ 

(3.104) 

From  Equations  (3.102)  and  (3.104),  it  can  be  easily  verified  that 

-  JC.^(-',)»il\Kb)  = - -  (3.105) 

irM  ^i„’(M) 

Hence  Equation  (3.103)  can  be  rewritten  as 


oo 


u°  = 


BllHkc/>')nl^(k^p) 


*‘’‘.=1  Bll\k„b)£;  \BP\k„b)-jC.^MB!!‘’(k^b) 


(3.106) 


l/=l/n 
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Due  to  the  attenuation  of  the  creeping  wave,  the  residue  series  solution  is  fast 
convergent  in  the  deep  shadow  region.  The  symbol  “oo”  is  used  to  emphasize 
that  only  a  few  dominant  modes  are  enough  to  obtain  the  accurate  result.  It  is 
found  that,  for  most  applications,  the  dominant  mode  in  Equation  (3.106)  is 
adequate  for  the  calculation  of  the  field  in  the  deep  shadow  region.  It  should  be 
noted  that  the  dominant  n^)^  mode  is  not  necessarily  the  first  mode  as  for  the 
conducting  cylinder.  A  detailed  discussion  on  the  dominance  of  the  residue  series 
mode  for  the  impedance  and  the  coated  cylinder  will  be  given  in  Chapter  IV. 

As  the  field  point  approaches  the  shadow  boundary  transition  region,  the  residue 
series  solution  becomes  slowly  convergent.  This  is  because  the  arc  length  of  the 
creeping  wave  propagation  is  short  in  the  transition  region  and  therefore  even  the 
modes  with  high  attenuation  can  no  longer  be  ignored. 

The  multiply  encircling  wave  can  also  be  obtained  from  the  residue  series 
solution  for  the  integral  representation  given  in  Equation  (3.15).  Thus, 

=  ±  v" _ _ 

J  i/=i/n 

e->n(2»+^)  +  g-ji'n(2ir-<l>) 

^  1  _  e-j2unr 

3.3.2  Diffracted  field 

Physically,  the  residue  series  representation  of  as  given  in  Equation  (3.106) 
can  be  interpreted  as  the  diffracted  field  which  explains  the  transmission  of  the 
wave  into  the  shadow  region.  However,  the  ray-picture  interpretation  can  not  be 
obtained  from  the  rigorous  residue  series  solutions  given  in  Equation  (3.106).  The 
Keller  type  GTD  ray  format  which  gives  the  residue  series  solution  a  creeping-wave 


(3.107) 


(a)  Impedance  Cylinder 


(b)  Coated  Cylind  er 
Figure  3.10:  Field  in  the  shadow  region 
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diffraction  interpretation  as  illustrated  in  Figure  3.10  can  be  obtained  by  following 
procedures.  From  the  Debye  approximation  of  the  Hankel  functions  are 


and 


where 


and 


Here,  we  make  assumptions  that  On  and  0n  approximated  by 


On  =  cos 


and 

Based  on  the  assumptions  given  in  Equations  (3.11^)  and  (3.113) 
kop'  sin  an  =  y/iKp')^  ~  ^  ^oy/p'^  -  ^  ==  ko»i 


(3.108) 


(3.109) 


(3.110) 


(3.111) 


(3.112) 


(3.113) 


(3.114) 


fcopsin/3„  =  yjikopf  -  «  ko\/p^-^  =  fco«2 


(3.115) 


As  shown  in  Figure  3.10,  aj  and  S2  are  the  distance  from  the  source  to  the  tan¬ 
gential  point  Qi  and  the  distance  from  the  tangential  point  Q2  to  the  field  point 
respectively.  Therefore, 


^Si\kop')  %  ^ I"-?  [*0^1  “  cos' 


(3.116) 


H^n{^op)  ^  I  -  t'n  cos' 


(3.117) 


Substituting  Equations  (3.116)  and  (3.117)  into  Equation  (3.106),  can  be  ex¬ 
pressed  as 


j  -2j  1  c-J*®*! 

St  — -  .  —  .  — . — —  ■  .  e  ■'2 

kob  y/koSi  y/koS2 


“oo” 

^  S  72)  a  (W  (2)  1  (3.118) 

n=i  H^JJikob)^  nr  {kob)  -  jC,^h{t')Hl  \kob)\ 

J  u=Vn 

As  shown  in  Figure  3.10,  9  is  defined  as 


0  =  <f>  —  cos' 


(3.119) 


Rearrangement  of  Equation  (3.118)  leads  to  Keller’s  surface  diffracted  ray  format 
for  2-dimensionsd  case. 


,  “00”  •  a  i  ^  t>~jko*2 

«  “'(Pi)  ■  E  •  e--’-""  ■  K'NO,)} 


(3.120) 
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(3.121) 


“Wi)-  E  {®4’*«?i)-«pI-(oUiAW  'p;'‘«?2)} 


where  a'^  and  represent  the  attenuation  constant  and  the  propagation  constant 
of  the  creeping  wave  respectively.  It  can  be  easily  shown  that 


(3.122) 


[y]  (3.123) 

^n^(Q)  is  referred  to  as  the  surface  diffraction  coefficients.  For  the  circular  cylin¬ 
der,  X>n^(Qi)  is  equal  to  Iln^*(<?2)*  From  a  comparision  of  Equation  (3.118)  and 
Equation  (3.120),  the  diffraction  coefficient  Hn^lQ)  can  be  expressed  as 


(3.124) 

Now  it  is  clear,  from  the  Equation  (3.120),  that  the  residue  series  solution  can 
be  interpreted  as  a  creeping-wave  diffraction  as  illustrated  in  the  Figure  3.10. 
It  should  be  noted  that  the  ray  path  associated  with  the  perfectly  conducting 
cylinder  remains  unchanged  for  both  the  impedance  and  the  coated  cylinder.  As 
shown  in  Figure  3.10,  the  incident  rays  which  strike  the  cylinder  at  the  tangential 
points((?i  and  Q2)  with  grazing  angle  launch  creeping  waves  which  travel  along 
the  arc  of  geodesic  path  on  the  cylinder  surface  following  Fermat  principle.  As 


these  creeping  waves  propagate,  they  decrease  exponentially  and  continuously  shed 
diffracted  rays. 

However,  in  using  the  ray  solution  given  in  Equation  (3.120),  one  has  to  be  aware  of 
the  approximations  given  by  the  Equation  (3.112)  and  (3.113).  For  the  conducting 
cylinder,  the  dominant  pole  is  determined  such  that  ss  kob.  Therefore,  the  above 
assumption  is  always  retisonable  if  the  source  or  the  field  point  is  not  very  close 
to  cylinder  surface,  where  GTD  ray  solution  fails.  However,  for  the  impedance  or 
the  coated  cylinder,  the  dominant  poles  are  not  always  close  to  Therefore, 
the  approximation  given  in  Equation  (3.112)  and  (3.113)  are  valid  only  when  both 
p’  and  p  are  large  compared  with  the  radius  of  the  cylinder  (i.e.  p'  >>  6  and 
p  >>  b).  Unless  both  p'  and  p  are  large,  a  non*ray  format  as  given  in  Equations 
(3.106)  should  be  used.  This  imposes  an  addtional  restriction  on  the  validity  of 
the  GTD  ray  solution  of  the  impedance  and  the  coated  cylinder. 

For  the  perfectly  conducting  cylinder,  the  expression  for  Vn^{Q)  can  be  further 
simplified  by  using  the  approximation  of  the  Hankel  function  in  terms  of  the  Miller 
type  Airy  function.  However,  for  the  impedance  and  the  coated  cylinder  where 
the  dominant  is  not  necessarily  close  to  kob,  Airy  functions  can  not  properly 
approximate  the  Hankel  functions  involved  in  the  expression  of  Vn^{Q). 

Ray  format  for  the  multiply  encircling  wave  can  also  be  obtaind  by  taking 
the  same  procedure  as  for  the  u®.  It  is  interesting  to  note  that 

®  ”  P=1  |B=-1 


Thus, 


-jkon 


»=1  n=l  V*o«2 
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-CO  “oo" 

+  «’(<?!)•  (3.126) 

Thus  consists  of  two  terms.  The  first  term  represents  multiple  encirclements 
in  the  counterclockwise  sense,  and  the  other  term  represents  those  in  the  opposite 
direction.  As  illustrated  in  Figure  3.11-(a),  u®  in  Equation  (3.106)  represents  the 
diffraction  field  which  travels  the  shortest  arc  length  on  the  cylinder  surface  in  a 
counterclockwise  direction  (i.e.  p  =  0).  Recall  that  ^  is  assumed  to  be  positive 
and  less  than  w. 
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(c)  p  =  1 

Figure  3.11;  Ray  encirclements  associated  with  the  summation  index  p 
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3.3.S  Equivalent  impedance  associated  with  the  creeping  wave  of  the 
coated  cylinder 

In  this  section  the  equivalent  impedance  associated  with  the  creeping  wave  of 
the  coated  cylinder  will  be  obtained. 

It  is  noted  C$  and  for  the  impedance  cylinder  are  given  as 

c,  =  |£  l«r  TM  ;  Cft  =  |i  for  TE  (3.127) 

Equivdently,  C,(t/n)  and  Ch(un)  for  the  coated  cylinder  may  be  expressed  as 

for  TM  ;  Ci(u„)  =  for  TE  (3.128) 

where  represents  the  equivalent  modal  surface  impedance  associated  with 

the  diffraction  field  of  the  coated  cylinder.  From  Equations  (3.128)  with  (3.8)  and 
(3.9)  we  obtain 

for  TMz 

(3.129) 

for  TEz 

(3.130) 

As  mentioned  in  the  previous  section,  the  equivalent  impedance  associated  with  the 
reflected  field  of  the  coated  cylinder  is  dependent  on  <^.  However,  these  equivalent 
modal  impedance j  associated  with  the  diffracted  held  as  given  by  Equation  (3.129) 
and  (3.130)  are  constant  along  the  cylinder  surface. 


ZlHvn)  =  j2, 


and 


ZlHon)  =  -jZi 


414l»)41’'(tia)  -  J?K4l‘)4;’'(tia) 
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3.4  Nuni«rical  results  and  discussion 
3.4.1  Numerical  results 

In  this  section,  the  GTD  solution  derived  in  this  Chapter  is  employed  to 
calculate  the  electromagnetic  field  scattered  from  the  impedance  and  the  coated 
cylinder.  The  geometry  of  the  problem  is  illustrated  in  Figure  1.1.  The  two- 
dimensional  line  source  is  located  at  =  0),  and  the  receiver  is  located  at 

F(p,^).  Plane  wave  incidence  is  the  case  when  p'  =  oo.  For  the  electric  line  source 
excitation  {TMg  case),  the  pattern  is  plotted  for  the  total  electric  field  intensity 
Eg,  and  for  the  magnetic  line  source  excitation  {TEg  case),  the  pattern  is  shown 
for  the  total  magnetic  field  intensity  Hg.  When  the  field  point  is  in  the  shadow 
region  as  shown  in  Figure  3.12-(a),  the  asymptotic  solutions  employ  two  ray  paths, 
one  is  OQ1Q2P  and  the  other  is  For  the  field  point  P  located  in  the  lit 

region  as  shown  in  Figure  3.12-(b),  the  ray  path  OQ1Q2P  is  replaced  by  the  path 
OQrP,  where  Qjp  is  the  specular  reflection  point  defined  by  $1  =  0^. 


Typical  results  for  the  field  patterns  due  to  a  line  source  radiating  in  the 

presence  of  a  impedance  cylinder  are  shown  in  Figure  3.13  —  3.16.  For  the  cylinders 

with  lossless  dielectric  coatings,  ttie  numerical  results  are  shown  in  Figure  3.17  ~ 

3.26.  The  field  patterns  due  to  a  plane-wave  incidence  on  the  cylinder  with  lossy 

dielectric  coatings  and  the  lossy  composite  coating  are  shown  in  Figure  3.27  ~ 
* 

3.36.  The  results  in  each  case  are  compared  with  eigenfunction  series  solutions. 
The  surface  impedances  of  the  impedance  cylinders  were  chosen  such  that  9  =  ±1 
where  q  is  defined  as 


q  =  -jmCg^h 


where 


(3.131) 


i.e. 
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(a)  Shadow  region 


Ray  paths  employed  for  the  GTD  solution 
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(3.132) 


q=-jm~  for  TM^ 

q=-jm^  for  TE^  (3.133) 

Zo 

Thus  the  surface  impedances  chosen  in  the  examples  are  pure  imaginary  (Positive 
q  represents  capcitive  impedance  for  TMz  and  inductive  impedance  for  TEz  and 
negative  q  represents  inductive  impedance  for  TM^  and  capacitive  impedance  for 
TEz.) 


3.4.2  Discussion 

The  numerical  results  of  GTD  solution  as  illustrated  in  Figure  3.13  ~  3.36 
show  excellent  agreement  with  the  eigenfunction  results  in  the  deep  lit  region  and 
in  the  deep  shadow  region.  However,  as  expected,  the  GTD  solution  fails  within 
the  shadow  boundary  (SB)  transition  region.  Neither  Geometrical  Optics  solution 
for  the  lit  region  nor  the  residue  series  solution  for  the  shadow  region  are  accurate 
in  the  vicinity  of  the  shadow  boundary.  These  two  solutions  have  a  discontinuity  at 
the  shadow  boundary.  The  fdlure  of  the  ordinary  GTD  solution  in  the  transition 
region  can  be  explained  as  follows.  In  the  formulation  of  the  G.O.  field  (i.e.,  the 
incident  field  u*  and  the  reflected  field  u*”),  an  assumption  is  made  that  the  Bessel 
functions  in  the  integrand  can  be  approximated  by  their  Debye  approximations. 
It  is  well  known  that  the  Debye  approximation  of  H},  ’  '(z)  is  appropriate  on  the 
condition  that 

\u  ~  z  \  >  \J\  (3.134) 
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As  the  held  point  approaches  the  shadow  boundary,  the  stationary  phase  point  i/« 
on  which  the  evaluation  of  the  integral  via  the  stationary  phase  method  is  based 
is  close  enough  to  the  argument  z  that  the  condition  in  Equation  (3.134)  is  not 
satisfied.  In  fact,  in  the  transition  repon  the  Bessel  functions  are  appropriately 
approximated  by  their  Watson’s  approximation  formula  where  the  Bessel  functions 
are  expressed  in  terms  of  the  Airy  functions.  Another  cause  for  the  failure  of  the 
ordinary  GTD  solution  in  the  transition  region  comes  from  the  fact  that  the  con¬ 
tribution  from  the  end  point  at  t/  =  koh  is  not  included  in  the  evaluation  of  the 
integral  via  the  stationary  phase  method.  In  the  transition  region  the  stationary 
phase  point  is  located  near  the  end  point  and,  therefore,  the  contribution  from 
the  end  point  can  no  longer  be  ignored.  The  residue  series  solution  which  is  fast 
convergent  in  the  deep  shadow  region  becomes  slowly  convergent  as  the  field  point 
approaches  the  shadow  boundary.  This  can  be  explained  by  the  attenuation  of  the 
creeping  wave  mode.  In  the  deep  shadow  region,  the  creeping  waves  attenuate  as 
they  travel  a  large  distance  in  terms  of  the  wave  length  along  the  surface  of  the 
cylinder.  Therefore  only  a  few  creeping  wave  modes  with  low  attenuation  constants 
contribute  to  the  field  in  the  deep  shadow  region.  However  in  the  transition  region, 
the  arc  distance  of  attenuation  of  the  creeping  wave  on  the  cylinder  suface  is  so 
small  that  even  the  creeping  wave  modes  with  high  attenuation  constant  cannot 
be  ignored.  Slow  convergence  of  the  creeping  wave  in  the  transition  region  makes 
it  difficult,  if  not  impossible,  to  obtain  the  result  with  acceptable  accuracy. 

The  failure  of  the  ordinary  GTD  solution  for  the  field  in  the  transition  region  ne¬ 
cessitates  the  development  of  Uniform  Geometrical  Theory  of  Diffraction  (UTD) 
which  will  be  discussed  in  the  Chapter  V, 
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•TO 


Figure  3.13:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 


Figure  3.15;  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 


•TO 


Figure  3.16:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 
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Figure  3,18:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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Figure  3.22;  Normalized  bist  atic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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Figure  3.23:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  ;  TMz  case 


Figure  3.24:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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Figure  3.25:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TMt  case 
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Figure  3.26:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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Figure  3.27:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  ;  TMz  case 
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Figure  3.28:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEt  case 
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Figure  3.29:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TMz  case 


Figure  3.30:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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Figure  3.31:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TM-  case 


Figure  3.32:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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Figure  3,33:  NormaJized  bistatic  scattering  pattern  of  coated  cylinder  ;  TMz  case 
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Figure  3.34:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz 


Figure  3.35:  NorniaJized  bistatic  scattering  pattern  of  coated  cylinder  :  TMz  case 
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Figure  3.36:  Normdized  bistatic  scattering  pattern  of  coated  cylinder  ;  TEz  case 
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CHAPTER  IV 


CHARACTERISTICS  OF  THE  CREEPING  WAVES 


In  the  previous  Chapter,  the  ray  format  for  the  diffracted  field  in  the  deep 
shadow  region  was  obtained.  It  is  well  known  that  the  geodesic  ray  path  of  the 
diffracted  field  includes  part  of  the  arc  along  the  cylinder  surface.  The  field  which 
propagates  along  the  arc  of  cylinder  surface  is  commonly  called  the  creeping  wave. 
It  is  the  purpose  of  this  Chapter  to  further  investigate  the  behavior  of  the  creeping 
wave  on  the  cylinder  surface.  As  mentioned,  the  effect  of  the  surface  impedance 
or  the  coating  of  dielectric/ferrite  material  is  more  obviously  illustrated  in  the 
diffracted  field  than  in  the  geometrical  optics  field  (incident  and  reflected  field). 
This  can  be  explained  by  the  fact  that  the  diffracted  field  is  launched  by  the  creep¬ 
ing  wave  which  propagates  on  the  cylinder  surface  under  direct  influence  of  the 
surface  impedance  or  the  coating  material.  Thus  it  is  important  to  investigate 
and  understand  the  fundamental  characteristics  of  the  creeping  wave  before  we 
attempt  to  control  the  diffraction  field  from  an  object  with  the  surface  impedance 
or  the  dielectric/ferrite  coating. 

For  the  perfectly  conducting  case,  the  behaviors  of  the  creeping  wave  on  the  curved 
surface  has  been  extensively  studied,  and  they  are  well  understood  now. 

For  the  impedance  cylinder,  this  problem  has  been  also  studied  by  Wait  [17], 
Streifer  [18]  and  Logan[20].  They  obtained  propagation  constants  for  the  creeping 
wave  based  on  the  Watson’s  approximation  where  the  Hankel  functions  are  ex¬ 
pressed  in  terms  of  the  Fock-type  Airy  functions.  In  this  Chapter,  their  results  are 
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reproduced  for  the  comparision  with  our  results.  Recently,  Wang  [4]  has  obtained 
more  accurate  numericed  data  for  the  propagation  constant  of  the  creeping  wave, 
the  Regge  poles,  and  the  natural  frequencies  of  an  impedance  cylinder  by  using 
both  the  Debye  and  the  Watson’s  approximations  of  the  Hankel  functions.  The 
creeping  wave  interpretation  for  the  resonance  and  the  correlation  between  the 
resonance,  the  Regge  poles,  and  the  natural  frequencies  of  the  impedance  cylinder 
has  been  also  discussed.  The  propagation  constants  for  the  first  creeping  wave 
mode  on  the  cylinder  with  lossless  coating  has  been  obtained  by  Elliott[19]. 

The  most  recent  work  on  this  problem  has  been  made  by  Paknys  and  Wang  [7,8]. 
They  made  an  important  contribution  by  giving  some  accuarate  numerical  results 
for  tliis  problem.  In  the  following  sections  of  this  Chapter,  some  part  of  their  work 
is  further  expanded  by  using  basically  the  same  numerical  approximations  for  the 
calculation  of  Hankel  functions  involved  in  the  formal  solution.  Some  aspects  per¬ 
taining  to  the  creeping  wave  which  were  not  treated  in  their  work  are  also  discussed 
in  this  Chapter. 

4.1  Propagation  and  Attenuation  Constants 

In  Chapter  III,  it  was  shown  that  the  ray  format  for  the  diffracted  field  is 
given  as 


=  «>(<?, )  ■  Pj'-Wi )  ■ 

=  Ti'((3,)  •»;•»((?,)  •  .  Vi\Q,)  (4.1) 

From  the  above  Equation  (4.1),  it  can  be  easily  seen  that 

=  ;  ^n  =  -R'[j]  (4.2) 
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Figure  4.1:  Creeping  wave 
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and  denote  the  attenuation  constant  and  the  propagation  constant  of  the 
creeping  wave  mode  respectively.  The  propagation  constant  determines  the  phase 
velocity  of  the  creeping  wave  and  the  attenuation  constant  represents  the  loss  of 
the  electromagnetic  energy  due  to  shedding  of  the  diffraction  field  and  the  loss  due 
to  the  surface  impedance  or  the  lossy  coating  material.  The  relative  phase  velocity 
of  the  creeping  wave  Vpjc  is  given  by 

c  Re[vn] 

where  c  denotes  the  velocity  of  the  light. 

The  propagation  constant  and  the  attenuation  constant  are  determined  from  the 
roots  of  the  transcendantal  equation  which  characterizes  the  impedance  or  the 
coated  cylinder  as  given  by 

=  0  (4.4) 

v=Vn 

The  complex  root  of  the  transcendental  equation  denoted  m  Vn  is  sometimes  called 
as  Reggt  pole.  It  is  well  known  that  there  exists  an  infinite  number  of  roots  for 
the  Equation  (4.4).  The  subscript  n  denotes  the  mode  number  of  the  root. 

4.1.1  Impedance  cylinder 

For  the  impedance  cylinder,  in  the  Equation  (4.4)  is  given  by 

Ct{u)  =  C,  =  Y^  for  TMj  (4.5) 

and 

Ch{o)  =  Ck  =  Y  for  TEz  (4.6) 
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Recent  work  by  Paknys  and  Wang  [7,8]  shows  a  trajectory  of  the  Regge  poles  for 

purely  imaginary  impedances.  In  their  work,  the  Hankel  functions  are  evaluated 

exactly  as  well  as  asymptotically  by  using  both  the  Debye  approximations  and 

the  Watson’s  approximation  according  to  the  ranges  of  validity.  The  accuracy 

of  their  results  is  fully  confirmed  by  the  excellent  agreement  between  the  field 

pattern  of  the  eigenfunction  series  solution  and  the  field  pattern  obtained  by  using 

their  propagation  constant.  Following  the  same  approximations  for  the  Hankel 

functions,  the  trajectories  of  the  Regge  poles  for  various  complex  impedances  are 

obtained  as  shown  in  Figure  4.2  for  kob=20.  Newton-Raphson’s  method  has  been 

used  for  the  numerical  search  for  the  roots  of  the  transcendental  equation. 

It  can  be  easily  seen  that  for  Cg  /,=0,  the  roots  of  transcendantal  equation  are 

(2)’ 

determined  from  the  well  known  zeros  of  HI  '  {koh)  which  are  commonly  called 
hard  poles  of  the  conducting  cylinder.  These  hard  poles  are  given  as 

( k  b\  IT 

i/n  =  Tjn  =K  kob  +  (  “I"  j  (4.7) 

where  /3„=1.0188,  3.2482,  4.8201, . 

For  the  small  value  of  the  hard  pole  of  the  conducting  cylinder  can  be 

used  as  initial  guessing  of  the  root  in  the  numerical  search  via  Newton-Raphson’s 
method. 

For  Cg^h  =  oo,  the  roots  of  transcendantal  equation  are  determined  from  the  zeros 
of  Hi,  {kob)  which  are  the  soft  poles  of  the  conducting  cylinder.  They  are  given 


where  an=2.3381,  4.0879,  5.5205, . 

Figure  4.2  illustrates  trajectories  of  Un  as  file  magnitude  of  C,  /,  increases  where 
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the  phase  of  the  Cg  f^  is  fixed.  The  trajectories  of  the  i^n  can  be  obtained  by  start¬ 
ing  from  the  hard  poles  rjn  and  slowly  increasing  the  magnitude  of 
When  the  phase  of  Cg  f^  is  between  —270®  ~  —90®,  the  impedance  has  a  negative 
real  value  and  therefore,  such  a  Cg  h  is  nonphysical.  Our  interests  are  the  cases 
where  the  phase  is  between  —90®  ~  90®. 

From  Figure  4.2,  we  can  observe  two  different  type  of  behavior  of  the  equiphase 
trajectories  for  the  first  mode  roots.  All  trajectoris  of  roots  of  the  first  mode  start 
from  the  first  hard  pole  which  is  denoted  as  t/j.  For  the  case  where  the  phase  of 
Cg  f^  is  between  70®  ~  90®  in  Figure  4.2,  the  real  part  of  first  mode  root  i/j  keeps 
increaing  as  the  magnitude  of  Cgj^  increases.  This  means  that  the  phase  velocity 
of  the  creeping  wave  becomes  slower  as  the  magnitude  of  C,  increases.  It  can  be 
also  noted  that  i/j  does  not  approach  the  first  soft  pole  We  define  such  a  root 
as  Elliott-type  root. 

For  the  case  where  the  phase  of  Cg^'^  lies  between  —90®  ~  65®  in  Figure  4.2,  the 
roots  approaches  the  soft  pole  denoted  as  as  the  magnitude  of  the  Cg^h  in¬ 
creases.  We  define  such  a  root  as  WaUon-type  root.  The  boundary  of  the  two 
different  types  is  somewhere  between  65®  and  70®. 

For  the  second  creeping  wave  mode,  all  the  roots  start  from  the  second  hard  pole 
r\2-  As  the  magnitude  of  Cg^h  increases,  the  equiphase  trajectories  exhibit  three 
different  types  of  behavior.  When  the  phase  of  C, /,  equals  90®  ~  70®,  the  root 
approaches  the  first  soft  pole  ^i.  When  the  phase  equals  65®,  the  real  part  of 
keep  increasing  as  the  maginitude  of  Cg^h  increases.  For  the  phase  of  Cg^}^  between 
60®  ~  —90®,  the  roots  approach  the  second  soft  pole  ^2-  Thus,  except  at  65®  where 
the  equiphase  trajectory  exhibits  Elliott-type  roots,  all  other  equiphase  trjectories 
of  the  second  mode  roots  belong  to  the  Watson-type. 

Wait  [17],  Streifer  [18]  and  Logan  [20]  also  exanhned  the  roots  of  the  tran- 
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g  = 

Un  from  Figure  4.2 

Un  from  Figure  4.3 

II 

1 

0 

21.919  -j0.38098 

22.121  -  j0.32760 

g=l.l  -  jU 

22.246  -  jO.28736 

22.496  -  jO.23492 

9=1.2- jO 

22.598  -  j0.20862 

22.905  -  j0.15897 

9=1.3- jO 

22.975  -  j0.14478 

23.355  -  jO.10006 

9=1.4- jO 

23.380  -  j0.09529 

23.849  -  j0.05769 

9=1.5  -  jO 

23.812  -  j0.05899 

24.390  -  j0.02987 

Table  4.1;  Cotnparion  of  roots 


Watson’s  approximation.  Sampled  data  of  i/j  for  9  =  1  —  jO  ~  9  =  1.5  -  jO 
obteuned  by  above  metod  are  compared  with  our  results  in  Table  4.1.  It  can  be 
observed  from  the  Table  4.1  that  there  are  significant  discrepancies  between  these 
two  results.  As  shown  in  Figure  4.4,  the  accuracy  of  these  propagation  constants 
are  also  tested  by  the  comparision  of  the  field  pattern  obtained  by  data  in  Table  1 
with  that  of  the  eigenfunction  solution  for  5  =  I  — jO*?®  1.2  — yO  and  5  =  1.3 -jO. 
It  can  be  observed  from  the  comparision  in  Figure  4.4  4.9  that  the  approximation 

of  the  Hankel  functions  as  given  in  Equation  (4.9)  is  not  appropriate  when  g  is 
larger  than  1.  This  is  due  to  the  fact  that  for  |»^n  ~  the  Debye 

approximation  of  the  Hankel  function  should  be  used  instead  of  the  approximation 
given  by  Equation  (4.9). 
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Figure  4.3:  The  roots  of  the  trascendentaJ  eouation  for  the  impedance  cylinder 

using  Equation  (4.12) 


98 


99 


- EIGENFUNCTION 


Figure  4.8:  Field  pattern  obtained  using  the  in  Figure  4.3 


6TD 


Figure  4.9:  Field  pattern  obtained  using  the  i/n  in  Figure  4.2 


4.1.2  Coated  cylinder 

For  the  coated  cylinder,  in  Equation  (4.4)  is  given  as 


C,{y)  =  -3 


.2,0 


Hl^\kib)Hl^\kia)  ~  Hl^\kib)Hi,'\kia) 


far  TMz 
(4.15) 


^  d^^\kib)Hl^Uha)  - 

^  Hl^\kib)dJ'^\kia)  -  Hl^\kib)Hl^'^' ikia) 


for  TEz 
(4.16) 


The  propagation  constants  of  the  creeping  wave  mode  for  the  cylinder  with  lossless 
dielectric  coating  was  first  obtained  by  Elliott  [19].  Recent  work  by  Paknys  and 
Wang  [7,8]  also  shows  the  trajectories  of  Un  for  the  changing  values  of  the  thickness 
of  the  lossless  dielectric  coating  and  the  accuracy  of  their  results  is  also  proved  by 
comparision  with  the  eigenfunction  results.  By  using  the  same  approximations  for 
the  Hankel  functions  as  employed  by  Paknys  and  Wang,  the  trajectories  of  for 
the  general  complex  tr  and  complex  Hr  are  obtained  and  shown  in  Figures  4.10  ~ 
4.13. 

From  the  Wronskian  of  the  Hankel  function,  it  can  be  easily  proved  that 


lim  Cs{y)  =  oo 

o— *6 


(4.17) 


lim  C/,(i/)  =  0 
a— *6 


(4.18) 


This  means  that  the  propagation  constant  for  the  coating  with  zero  thickness  is 
the  same  as  for  the  perfectly  conducting  cylinder  and  therfore,  the  well  known 


soft  and  hard  poles  of  conducting  cylinder  can  be  used  as  the  initial  guess  of 
roots  for  the  thin  coating  in  the  numerical  search  of  wn  via  Newton-Raphson’s 
method.  From  numerical  results  given  in  Figure  4.10  ~  4.13,  we  can  also  observe 
two  different  types  of  behavior  in  the  trajectories  as  the  thickness  of  the  coating 
increases.  When  the  real  part  of  Vn  increases  (i.e.  the  phase  velocity  decreases)  as 
the  coating  thickness  increases,  we  refer  to  this  trajectory  as  Elliott- type.  It  can  be 
noted  that  the  Elliott-type  trajectory  is  obtained  only  when  the  coating  material 
has  a  low  loss  tangent.  For  the  material  coating  with  a  high  loss  tangent,  the 
trajectory  of  i^n  represents  the  Watson-type.  It  can  be  observed  that  the  Watson- 
type  i/n  approaches  a  limiting  point  for  each  creeping  wave  mode  as  the  thickness 
of  coating  increases.  If  the  coating  material  has  a  high  loss  tangent,  the  field  inside 
the  coating  material  decreases  exponentially  as  the  field  point  approaches  the  inner 
conducting  cylinder.  Thus  the  field  inside  the  coating  material  is  confined  to  near 
the  surface  of  the  cylinder.  For  this  case,  further  increase  of  the  coating  thickness 
does  not  affect  the  propagation  and  the  attenuation  constant  of  the  creeping  wave 
mode.  Figure  4.14  shows  the  relative  phase  velocity  and  the  attenuation  of  the 
first  creeping  wave  mode  versus  coating  thickness.  Curves  are  shown  for  four 
different  sizes  of  the  coated  cylinder  {kob=20,  40  ,100,  oo)  with  coating  material  of 
=  5.1513  —  j4.253  and  pr  =  1-  The  attenuation  shown  in  Figure  4.14  represents 
the  attuation  per  meter  at  a  frequency  of  10  GHz.  The  creeping  wave  on  the 
cylinder  with  infinite  radius  (i.e.  kob  =  oo)  is,  in  fact,  the  surface  wave  on  the 
planar  ground  slab^  with  the  same  coating  material. 


’■The  numerical  data  for  the  surface  wave  were  obtained  by  Prof.  Richmond  of  the  Ohio  State 
University. 
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Figure  4.10:  The  roots  of  the  transcendental  equation 
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Figure  4.11:  The  roots  of  the  transrendental  equation 
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Figure  4.12:  The  roots  of  the  transcendental  equation 
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Figure  4.13:  The  roots  of  the  transcendental  equation 
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4.2  Dominance  of  Creeping  Wave  Modes 

It  is  known  that  the  coated  cylinder  supports  an  infinite  number  of  creeping 
wave  modes.  In  this  section,  the  contribution  of  the  different  creeping  wave  inodes 
to  the  diffraction  field  is  compared.  The  residue  series  solution  for  the  diffraction 
field,  which  is  obtained  from  the  Watson’s  transformation  of  slowly  convergent 
angular  eigenfunction  series  solution,  is  fast  convergent  in  the  deep  shadow  region. 
Thus  in  the  actual  calculation  of  the  diffracted  field,  one  or  two  dominant  modes 
are  enough  to  obtain  a  numerically  accurate  result.  It  is  important  therefore  to 
determine  which  mode  is  dominant  over  other  modes.  One  should  not  conclude 
that  a  certain  mode  is  dominant  over  other  modes  simply  because  that  it  has  a 
smaller  attenuation  constant.  This  is  because  the  magnitude  of  the  diffracted  field 
is  dependent  not  only  on  the  attenuation  constant  but  also  on  the  diffraction  co¬ 
efficient  given  in  Equation  (3.124).  For  the  perfectly  conducting  cylinder,  the  first 
creeping  wave  mode  is  always  dominant  over  all  other  modes  and  therefore  only 
the  first  mode  is  adequate  in  practice.  This  is  partly  because  the  first  mode  has  a 
lower  attenuation  constant  than  other  higher  modes.  However,  this  is  not  always 
true  for  the  impedance  or  the  coated  cylinder.  For  the  coated  cylinder,  the  domi¬ 
nance  of  the  creeping  wave  mode  changes  as  the  thickness  of  the  coating  changes. 
Thus,  for  some  thickness  of  coating,  the  second  or  even  higher  mode  has  a  smaller 
attenuation  constant  and  a  larger  diffraction  coefficient  than  the  first  mode  and 
therefore  is  dominant  over  the  first  mode.  This  fact  is  well  demonstrated  by  fol¬ 
lowing  example. 

Figure  4.15  shows  the  nomalized  backscattering  width  of  the  coated  cylinder.  The 
radius  of  the  inner  conducting  cylinder  is  35/27r  Xo  (i.e.,  ^<,6=35).  A  lossless  dielec¬ 
tric  material  with  €r=2.56  and  pr=l  is  coated  on  the  conducting  cylinder.  The 
incident  plane  wave  has  a  magnetic  field  parallel  to  the  cylinder  axis  (i.e.  TEz 
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case).  It  is  known  that  the  distinctive  ripples  in  the  backscattered  field  is  due  to 
the  creeping  wave  resonance  [4].  Furthermore,  the  resonances  in  Figure  4.15  can 
be  traced  to  the  contribution  of  two  creeping  wave  modes. 

For  a  very  thin  coating,  the  diffraction  field  is  small  due  to  high  attenuation  and 
therefore  normalized  backscattering  width  is  mainly  determined  by  the  reflection 
field.  As  the  coating  material  is  lossless,  normalized  backscattering  width  is  zero 
dB  for  very  thin  coating.  The  resonance  comes  from  the  first  creeping  wave  mode 
whose  trajectory  of  propagation  constant  (i.e.  i^i)  as  a  function  of  coating  thick¬ 
ness  is  shown  in  Figure  4.17.  The  resonance  from  the  first  creeping  wave  mode 
becomes  negligible  after  the  thickness  of  the  coating  exceeds  0.2  Xq-  For  coating 
thicknesses  of  0.2  Aq  ~  0.45  A©,  creeping  wave  is  very  weak  compared  with  the 
reflected  field  and  thus  has  a  negligible  contribution  to  the  total  backscattering 
field.  As  the  coating  thicknesss  increases  further,  the  resonance  originated  from 
the  second  creeping  wave  begins  to  appear  in  the  backscattering  width  between 
d  =  O.45A0  and  d  —  0.7Ao.  The  trajectory  of  the  propagation  constants  for  the 
second  creeping  wave  mode  (i.e.  1/2)  is  also  plotted  in  Figure  4.17.  Thus,  the  sec¬ 
ond  creeping  wave  mode  plays  an  important  role  in  the  total  backscattering  field 
for  coating  thickness  of  0.45  Aq  ~  0.7  A©  and  disappears  for  further  increase  of 
the  coating  thickness.  It  is  also  expected  that  if  we  further  increase  the  coating 
thickness,  the  third  or  even  higher  modes  can  be  dominant.  This  phenomena  can 
be  explaned  as  follows. 

The  magnitude  of  the  diffracted  field  depends  on  two  factors;  the  diffraction  coef¬ 
ficient  and  the  attenuation  constant  of  the  creeping  waves.  It  turns  out  that,  for 
the  range  from  d  =  0.2Ao  to  d  =  0.45Ao,  even  the  attenuation  is  small  (see  Figure 
4.17),  the  diffraction  coefficient  for  the  first  creeping  wave  mode  as  shown  in  Figure 
4.16  is  insignificant  to  contribute  to  the  backscattered  field.  In  the  same  range, 
the  contribution  to  the  backscatered  field  due  to  the  second  creeping  wave  mode  is 
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also  minimal  due  to  the  attenuation.  However,  when  the  thickness  of  the  coating  is 
increased  beyond  0.45Ao,  the  attenuation  rate  for  the  second  creeping  wave  mode 
drops  significantly  such  that  its  contribution  dominates  the  backscattered  field.  It 
is  of  interest  to  point  out  that  the  thickness  of  the  dielectric  coating  for  the  coated 
cylinder  at  which  the  resonances  begin  to  occur  corresponds  approximately  to  the 
cutoff  frequencies  of  the  surface  waves  supported  by  the  planar  grounded  slab  {7j. 
Similar  creeping  wave  resonances  occur  in  the  backscattering  width  for  the  TMz 
case  as  shown  in  the  Figure  4.18.  The  corresponding  trajectories  of  the  propaga¬ 
tion  constants  are  shown  in  the  Figure  4.20  and  the  diffraction  coefficients  for  the 
TMz  modes  are  shown  in  Figure  4.19.  The  explanation  for  the  TEz  resonances 
also  applies  to  the  TMz  case.  However,  it  is  found  that  the  second  resonance  oc- 
curing  between  d  =  0.6Ao  and  d  =  0.75Ao  as  shown  in  Figure  4.18  comes  from  a 
creeping  wave  with  a  different  origin.  Refering  to  Figure  4.20,  it  is  seen  that  as  the 
thickness  of  the  coating  is  decreased,  the  trajectory  of  this  unusual  mode  does  not 
lead  to  the  usual  second  creeping  wave  mode  for  the  perfectly  conducting  cylinder 
(i.e,  ^2)-  appears  that  there  exist  an  additional  set  of  creeping-wave  modes 
for  cylinders  with  dielectric  coating.  These  ‘new’  creeping-wave  modes  yield  the 
dominant  contribution  to  the  diffracted  field  when  the  coating  becomes  thick. 


Ill 


DIFFRACTION  COEFFICIENT  IDI 


Figure  4.17;  Propagation  and  the  attenuation  constant 


113 


DIFFRACTION  COEFFICIENT  IDI  r /irb 


k,*  •  ss 


-4-1 


Figxire  4.18:  Normalized  echo  width 
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4.3  Creeping  waves  on  a  large  cylinder 

In  this  section,  we  investigate  the  beljavior  of  creeping  waves  on  the  coated 
surface  with  electrically  large  radius  of  curvature.  The  investigation  in  this  section 
is  prompted  by  an  apparent  discrepancy  between  the  attenuation  constants  for  the 
creeping  waves  on  a  coated  curved  surface  and  those  for  the  surface  wave  on  a  pla¬ 
nar  grounded  slab.  Referring  to  Figure  4.14,  as  the  circumference  of  the  circular 
cylinder  is  increased  from  20  to  100,  the  attenuation  constants  for  the  creeping 
waves  do  not  approach  the  planar-slab  limit.  This  contradicts  past  experience  for 
the  perfectly  conducting  cylinder  with  a  lossless  coating.  The  radius  of  the  circu¬ 
lar  cylinder  is  increased  to  a  even  larger  value  in  order  to  check  the  limiting  case. 
The  attenuation  constants  were  calculated  for  several  cylinders  with  large  radius, 
and  interesting  results  were  obtained.  Figure  4.21  presents  the  the  roots  of  the 
transcendental  equation  of  the  circular  cylinder  [b  =  50Ao,  i.e.  kob  =  IOOtt)  with 
a  high  loss  material  coating  (cr  =  5.1513  -  j4.253,  pr  =  1  ~  jO).  In  Figure  4.21, 
the  trajectories  Un  for  the  first  five  of  the  infinite  creeping  wave  modes  are  shown 
as  a  function  of  the  thickness  of  coating  from  zero  to  a  final  value  of  0.1  Ac.  It 
is  observed  from  Figure  4.21,  that  the  imaginary  part  of  i/n  of  the  first  creeping 
wave  mode  becomes  larger  than  those  of  mode  2,  3,  and  4  as  the  thickness  of  the 
coating  is  increased.  The  corresponding  attenuation  constants  are  shown  in  Fig¬ 
ure  4.22.  It  can  be  seen  that  the  attenuation  constant  for  the  first  creeping  wave 
mode  behaves  like  that  of  the  surface  wave  associated  with  the  planar  grounded 
slab  as  shown  in  Figure  4.14.  The  fields  associated  with  the  first  mode  are  rapidly 
attenuated  to  a  negligible  value.  This  phenomena  does  not  occur  for  the  smaller 
cylinder  (kob  =  100  i.e.  b  =  l00/2w  Xg  Figure  4.14).  The  diffraction  coefficient 
of  each  mode  is  also  plotted  in  Figure  4.23.  The  dominance  of  the  mode  can  be 
determined  by  comparing  both  the  attenuation  and  the  diffraction  coefficient  of 
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each  mode. 

To  further  demonstrate  the  limiting  case,  numerical  results  obtained  for  a  cylinder 
with  b=100  \o  are  shown  in  Figure  4.24  and  4.25.  It  is  seen  that  the  attenua¬ 
tion  constants  of  the  first  creeping  wave  mode  for  both  cylinders  (b=50  Xo  and 
b=100  Xq)  converge  to  each  other  and  both  approach  the  planar  limit.  Finallly, 
the  diffraction  coefficients  for  the  creeping  wave  modes  of  the  coated  curved  sur¬ 
face  with  a  radius  of  100  Xo  are  shown  in  Figure  4.26.  Thus,  in  this  section  we 
demonstrated  the  fact  that  the  first  creeping  wave  mode  behaves  the  same  way  as 
the  surface  wave  associated  with  the  planar  grounded  slab. 
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Figure  4.21;  Trajectory  of  roots 
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Figure  4.22:  Attenuation  per  meter 
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Figure  4.24:  Trajectory  of  roots 
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Figure  4.26;  Diffraction  coefficient 


CHAPTER  V 


FIELD  IN  THE  TRANSITION  REGION 


As  seen  from  the  numeriacl  results  of  examples  given  in  Chapter  III,  the  GTD 
solution  fails  within  the  shadow  boundary  (S.B.)  transition  region.  The  various 
causes  for  this  failure  are  also  explained  in  Chapter  III.  In  this  Chapter,  we  at¬ 
tempt  to  correct  the  failure  of  the  ordinary  GTD  in  the  transition  regions. 

For  a  perfectly  conducting  cylinder,  the  uniform  GTD  solution  which  is  often  ab¬ 
breviated  as  UTD  has  already  been  developed  and  successfully  applied  to  many 
engineering  problems.  UTD  solution  for  the  curved  surface  was  first  introduced  by 
Pathak  [3]  for  the  canonical  problem  of  plane  wave  scattering  by  a  smooth  perfectly 
conducting  circular  cylinder.  As  its  name  implies,  uniform  GTD  is  “uniform”  in 
the  sense  that  it  remains  valid  within  the  transition  regions  and  it  automatically 
reduces  to  the  ordinary  GTD  exterior  to  the  transition  regions  where  the  ordinary 
GTD  becomes  valid  and  accurate.  Both  GTD  and  UTD  employ  the  same  ray 
pathes  and  thus  both  solutions  have  the  same  ray  formats. 

For  the  impedance  cylinder,  UTD  solution  has  already  been  developed  by  a  heuris¬ 
tic  extension  of  the  UTD  solution  for  the  conducting  cylinder  [5].  For  the  cylinder 
with  dielectric/ferrite  coating,  the  solution  which  gives  a  uniform  result  in  the 
transition  region  has  not  yet  been  developed.  However,  the  fact  that  both  the 
impedance  and  the  coated  cylinder  employ  the  same  ray  format  outside  the  tran¬ 
sition  region  as  for  the  conducting  cylinder  suggests  the  possibility  that  the  same 
uniform  GTD  solution  developed  for  the  conducting  cylinder  might  also  be  used 
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for  the  coated  cylinder  within  the  transition  region.  Based  on  this  conjecture,  we 
heuristicadly  generalize  the  UTD  solution  for  the  perfectly  conducting  cylinder  to 
both  the  impedance  and  the  coated  cylinder  with  necessary  modifications  in  the 
transition  function  due  to  the  differences  in  the  boundary  conditions.  The  diffici- 
culty  in  finding  the  uniform  solution  for  the  impedance  and  the  coated  cylinder  is 
that  the  numerical  data  for  the  transition  function  F, which,  is  well  tabulated 
for  the  conducting  cylinder  is  not  available.  Instead  of  direct  numerical  evaluation 
of  the  transition  function,  an  exact  eigenfunction  solution  together  with  the  GTD 
solution  is  used  to  obtain  the  necessary  numerical  data  for  The  validity  and 

limitation  of  the  heuristic  method  has  also  been  tested  by  some  examples. 
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5.1  Uniform  GTD  solution  for  the  impedance  and  the  coated  cylinder 


In  this  section  we  generalize  the  uniform  GTD  solution  for  the  conducting 
cylinder  to  the  impedance  and  the  coated  cylinder  based  on  the  conjecture  that 
the  ray  format  given  for  UTD  solution  of  the  conducting  cylinder  can  be  retained  to 
both  the  impedance  and  the  coated  cylinder.  The  lengthy  asymptotic  analysis  and 
the  ray  format  of  UTD  solution  for  the  conducting  cylinder  are  given  by  Pathak 
[3]. 

5.1.1  Lit  region 


From  the  result  given  for  the  conducting  cylinder,  the  ray  format  for  the 
impedance  and  the  coated  cylinder  in  the  lit  region  can  be  obtained  by  replacing 
P{(']  with  P{(',q).  Thus  when  the  field  point  Pi  is  in  the  lit  region,  u°  is  given 
as 


+  h 


(5.1) 


where  the  reflected  ray  caustic  distance  />*■  can  be  calculated  from  the  relationship 


i  =  f  +  (5.2) 

p*"  6cos^j 

The  subscript  s  and  h  denotes  TMz  (soft)  and  TEz  (hard)  case  respectively.  In 
Equation  (5.1),  the  surface  reflection  coefficient  72,  in  the  ordinary  GTD  solution 
is  now  replaced  by  the  generalized  reflection  coefficient  defined  as 


(■5.3) 


where 
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Figure  5.1:  Field  in  the  Hi  part  of  transition  region 

(k 

i'  =  -2m  cos ;  m  =  I  ^  j  (5.4) 

and 

=  2ko  cos2^,  (5.5) 

The  relevant  parameters  $i,  /i ,  I2  and  the  specular  reflection  point  Qji  are  shown 
in  Figure  5.1. 

The  transition  function  F(t)  which  involves  a  Fresnel  integral  is  defined  as 


127 


(6.6) 


F(i)  =  2j-\/x  e-^*  /y_  c  dr 
Jy/z 

The  transition  function  F(x)  plays  a  key  role  in  enforcing  the  continuity  and 
boundedness  of  the  total  field  at  the  shadow  boundary. 

The  transition  function  which  is  generally  referred  as  the  Pekeris  cant 

function  is  defined  as 


V'jr.)  -  qV{r) 
J-oo-j^Wl^(T)-qW2{T) 


(5.7) 


where  q  is  defined  as 


q  =  -jmC^^hit')  (5-8) 

Thus  for  the  impedance  cylinder,  ^  is  a  constant  given  as 

q  =  -jm^  for  TMz  :  9  = -jm  ^  for  TEz  (5.9) 

However  for  the  coated  cylinder,  9  is  a  function  of  u.  For  the  simplicity  of  the 
calculation,  we  assume  that  q  can  be  approximated  as 


q  %  -jmC,,fc(i/„)  (5.10) 

where  !/„  is  the  dominant  n^/,  root  of  the  transcendantal  equation  given  in  Equation 
(4.4).  The  Fock-type  Airy  functions  V’(r)  and  1V2(t)  are  defined  as 

2jT(t)  =  H^’(r)- W^2(t)  (5.11) 


and 
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(5.12) 


2  =  4=  /  dr 


The  contour  Fj  runs  from  oo  to  oo  —  je,  and  r2  is  the  complex  conjugate 

of  Fj.  The  transition  function  can  eJso  expressed  in  terms  of  the  the 

Fock  integrals  p*{x,q)  and  g*(x,q). 


1 

2xy/ir 


(5.13) 


5.1.2  Shadow  region 

A  ray  format  for  the  shadow  region  can  also  be  obtained  from  the  result  given 
for  the  conducting  cylinder  by  replacing  P{()  with  P((,q).  When  the  field  point 
is  in  the  shadow  region,  the  ray  format  for  u®  is  given  by 

p-jkoS2 

n^{Pd)^AQl)Ts,h—j^  (5.14) 

From  Equation  (3.121),  the  transfer  function  T, /,  for  the  ordinary  GTD  solution 
is  given  as 


Kh  S  i-  E  )  «P[-(»'n  +  0nWK'‘iQ2)  (5.15) 

v^o  n 

where  the  diffraction  coefficient  Vn^iQ)  is  given  by  Equation  (3.124).  However,  in 
the  unform  GTD,  the  transfer  function  Tg  f^  is  replaced  with  a  generalized  transfer 
function  7^  ^  defined  as 


Psfh 


~m  € 


-jkoi 


e-J"’/'*  r  rf  1  * 


(5.16) 
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Figure  5.2;  Field  in  the  shadow  part  of  transition  region 
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where 


^  =  ma  ;  m 


(5.17) 


and 


Si  +  S2  2m^ 


(5.18) 


i  =  b  a 


(5.19) 


The  relevant  parameters  si,  S2,  t,  ol  and  two  tangent  points  Qi,  Q2  are  shown 
in  Figure  5.2.  The  uniform  GTD  solution  reduces  to  the  ordinary  GTD  solution 
in  the  deep  shadow  region. 
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5.2  Transition  function 


The  transition  function  has  also  been  studied  by  Wait  and  Conda 

[17]  for  a  cylinder  with  finite  conductivity.  The  quantity  >6,  in  fact, 

equivalent  to  the  G(^)  defined  in  [17,  eq(22)].  For  a  cylinder  with  finite  conduc¬ 
tivity,  the  function  G(^)  has  been  evaluated  via  direct  numerical  integration  and 
tabulated  for  the  various  values  of  g.  However,  for  a  positive  real  value  which 
occurs  for  a  cylinder  with  a  pure  reactive  surface  impedance  or  with  a  lossless 
coating  no  numerical  data  have  been  given.  Under  this  condition, the  direct  nu- 
mericad  integration  is  more  involved  because  the  integrand  has  a  singularity  close 
to  the  integration  path.  Therefore  a  heuristic  approach  [5],  rather  than  the  direct 
integration,  is  employed  in  this  section  to  deduce  numerical  results  for  the  transi¬ 
tion  integral.  This  engineering  approach  is  described  in  following  paragraph. 


a)  Lit  part  of  the  transition  region 
From  Equation  (5.1)  and  (5.3)  we  obtain 


J 


>  - 


F(X^) 


€. . e-M 

P”  +  h 


Note  that 


(5.20) 


(5.21) 


Therefore,  the  Fock  integrsJ  can  be  expressed  as 


J 


4.___ir(A-i) 
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(5.22) 

For  the  numerical  value  of  u*,  we  use  the  eigenfunction  solution  given  in  Chapter 
I  and  for  we  use  the  residue  series  solution  given  in  Equation  (3.12C). 

a)  Shadow  part  of  the  transition  region 

For  the  shadow  part  of  the  transition  region,  the  same  procedure  employed  for  the 
lit  part  can  be  used  to  obtain  the  numerical  data ’for  the  From  Equation 

(5.14)  and  (5.16) 

=  u’(Pj)  -  (5.24) 


■  K  e-j^o»2 

(5.23) 


Thus,  we  obtain 


uHPd)  -  ^^^(Pd) 


e->V4 


(5.25) 


As  before,  we  calculate  u^{P^)  from  the  exact  eigenfunction  series  solution  and 
from  the  residue  series  solution. 


For  the  impedace  cylinder,  the  numercal  data  for  the  Fock  integral  9) 

or  obtained  from  the  Equation  (5.22)  and  (5.25)  are  given  in  Figure 

5.3  and  5.4  for  several  different  q.  It  can  be  observed  from  Figure  5.3  and  5.4  that 
the  numerical  results  of  the  Fock  integral  are  weakly  dependent  on  kob,  kop'  and 
kop.  Therefore,  curves  shown  in  Figure  5.3  and  5.4  can  be  used  as  “universal” 
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curves,  to  calculate  the  scattered  field  of  the  impedance  cylinder  with  arbitrary 
radius. 

For  the  coated  cylinder,  the  numercal  data  for  the  Fock  integral  obtained 
using  the  Equation  (5.22)  and  (5.25)  are  plotted  in  Figures  5.5,  5.6,  5.8  and  5.9 
for  various  combination  of  ei  and  kob.  The  Fock  integrals  obtained  for  the  TMx 
coated  cylinders  sire  shown  in  Figures  5.5  and  5.5  and  those  for  the  TEi  coated 
cylinders  are  shown  in  Figures  5.8  and  5.9.  Note  that  the  Fock  integral  depends 
on  the  q-parameter.  The  q-parameter  which  is  determined  only  by  the  constant 
surface  impedance  Zg  for  the  impedance  cylinder  is  now  a  function  of  the  radius 
of  coated  cylinder,  the  thickness  and  the  material  of  the  dielectric  coating  for  the 
coated  cylinder.  However,  it  is  found  that,  for  the  thin  coating  case,  the  transition 
integrals  obtained  are  almost  identical  providing  the  q-values  are  the  same.  This 
can  be  seen  from  the  results  presented  in  Figures  Figures  5.5,  5.6,  5.8  and  5.9. 
It  turns  out  that,  providing  both  have  the  same  q-value,  the  transition  integrals 
obtained  for  the  coated  circular  cylinder  compare  surprisingly  well  with  those  ob¬ 
tained  for  the  impedance  cylinder.  This  fact  can  be  established  by  comparing  the 
results  shown  in  Figures  5.7  and  5.10  with  those  shown  in  Figures  5.5,  5.6  and 
5.8,  5.9  respectively.  As  a  consequence,  the  transition  integral  obtained  for  the 
impedance  cylinder  could  be  employed  in  the  uniform  GTD  solution  to  calculate 
the  scattered  fields  in  the  transition  region  for  the  coated  cylinder.  The  advan¬ 
tage  here  is  that  the  procedure  involved  to  deduce  the  transition  integral  for  the 
impedance  cylinder  is  much  less  complicated  than  that  employed  for  the  coated 
cylinder.  Note  that  the  above  comments  only  apply  to  coated  cylinder  with  a  thin 
coating.  At  this  stage,  it  is  not  yet  clear  whether  the  concept  of  uniform  GTD 
solution  with  ‘universal'  transition  integral  is  applicable  for  the  case  of  a  cylinder 
covered  with  a  thick  dielectric  coating.  This  needs  to  be  further  investigated. 
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Figure  5.5:  The  transition  integral  :  Coated  Cylinder 
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Figure  5.6:  The  transition  integral  e  :  Coated  Cylinder 
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Figure  5.8;  The  transition  integral  ;  Coated  Cylinder 
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Figure  5.9:  The  transition  integraJ  e  :  Coated  Cylinder 
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Figure  5.10:  The  transition  integral  e  or  ;  Impedance 

cylinder 


5.3  Numerical  results  and  discussion 

In  this  section,  the  UTD  solution  is  employed  to  calculate  the  electromagnetic 
field  scattered  from  the  impedance  and  the  coated  cylinder.  The  geometry  of 
the  problem  and  the  ray  paths  are  the  same  as  those  of  GTD  solution.  It  is 
noted  that  the  results  obtained  with  the  GTD  solution  given  in  Chapter  III  show 
excellent  agreement  with  the  eigenfunction  results  in  the  deep  lit  and  the  deep 
shadow  regions.  However  the  GTD  solution  fails  in  the  transition  region  adjacent 
to  the  shadow  boundaries.  The  angular  range  of  the  regions  where  the  pure  ray- 
optical  solution  fails  depends  on  the  range  where  the  approximations  involved  in 
the  formulation  of  ordinary  GTD  solution  fails.  From  the  numerical  results  given 
in  Chapter  III,  It  can  be  observed  that  the  lit  part  of  the  transition  region  has 
wider  angular  range  than  the  shadow  part. 

It  is  well  known  that  the  uniform  GTD  reduces  to  the  ordinary  GTD  in  the  deep 
lit  and  the  deep  shadow  region.  Thus,  in  obtaining  the  numerical  results  given  in 
Figure  5.11  ~  5.34,  the  uniform  GTD  solutions  givexr  in  Equation  (5.1)  and  (5.14) 
are  employed  only  for  the  transition  region.  For  the  deep  lit  and  the  deep  shadow 
region  the  ordinary  GTD  solutions  given  in  Chapter  III  are  used.  The  examples 
given  in  Figure  5.11  ~  5.34  are  chosen  to  have  the  same  configurations  as  those 
given  in  Figure  3.13  ~  3.36  for  the  comparision  of  the  GTD  and  the  UTD.  It  can 
be  observed  that  the  uniform  GTD  solution  provides  a  continuous  pattern  across 
the  shadow  boundaries  and  compares  well  with  the  eigenfunction-solution  results 
in  the  transition  region.  Due  to  the  numerical  difficulty  in  locating  the  high  order 
dominant  root  of  the  tracendental  equation,  the  examples  given  in  Figure  5.11 
~  5.34  are  confined  to  thin  material  coatings  where  the  dominant  creeping  wave 
mode  is  either  the  first  or  the  second  mode.  At  this  stage  it  is  not  clear  w'hether 
the  solution  is  applicable  to  the  thick  coating  case. 


Figure  5.11:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 


Figure  5.12:  Normalized  bistatic  scattering  pattern  of  impedance  cylinder 
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Figure  5.17:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TM2  case 


Figure  5.18:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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Figure  5.21:  Normalized  bistaiic  scattering  pattern  of  coated  cylinder  :  TMz  case 


Figure  5.22:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TE:  case 
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Figure  5.25:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TMz  case 


Figure  5.27:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TM-  case 


Figure  5.28:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 


Figure  5.30:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TE-  case 
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Figure  5.31:  Normalized  bistatic  scatteriiig  pattern  of  coated  cylinder  :  TM-  case 


154 


Figure  5.33:  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TM^  case 


Figure  5.34;  Normalized  bistatic  scattering  pattern  of  coated  cylinder  :  TEz  case 
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CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 


An  approximate  asymptotic  high-frequency  solution  is  obtained  for  the  field 
exterior  to  a  two-dimensional  circular  conducting  cylinder  with  a  thin  dielec¬ 
tric/ferrite  coating.  The  angular  eigenfunction  solution  of  the  problem  is  obtained 
from  the  two-dimensional  line-source  Green’s  function  subject  to  the  boundary 
conditions  on  the  cylinder  surface  and  the  Sommerfeld  radiation  condition  at  the 
infinity.  The  eigenfunction  solution  is  slowly  convergent  when  the  radius  of  the 
cylinder  is  large  in  term  of  the  wave  length.  Furthermore,  the  physical  interpreta¬ 
tion  of  the  scattering  machanism  is  not  possible  from  the  eigenfunction  solution. 
In  order  to  obtain  the  ray  optical  solution  of  the  problem,  the  angular  eigenfunc¬ 
tion  is  transformed  into  the  integral  representation  and  the  integral  is  evaluated 
asymptotically  via  the  stationary  phase  method  for  the  field  in  lit  region  and  via 
the  residue  theorem  for  the  field  in  shadow  region.  The  obtained  solution  is  cast 
in  the  format  of  the  geometrical  theory  of  diffraction  (GTD). 

In  the  lit  region,  the  geometrical  optics  (GO)  solution  consists  of  the  direct  inci¬ 
dent  ray  and  the  reflected  ray.  The  ray  picture  sissociated  with  the  reflected  ray  for 
the  coated  cylinder  remains  the  same  as  that  for  the  conducting  cylinder  without 
coating.  The  only  difference  is  that,  for  the  coated  surface,  the  specular  reflection 
point  is  located  on  the  coating-air  interface.  It  is  shown  that  for  a  thin  coating 
case,  the  reflection  coefficient  for  the  field  reflected  at  the  specular  reflection  point 
on  the  coated  curved  surface  can  be  approximated  by  the  reflection  coefficient  for 
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a  grounded  planar  dielectric/ferrite  slab  with  a  plane  wave  incidence. 

In  the  shadow  region,  the  residue  series  solution  can  be  interpreted  as  the  usual 
creeping  wave  diffraction  of  the  GTD  solution  with  the  same  ray  path  as  the 
diffracted  field  associated  with  the  conducting  cylinder  without  coating.  The  prop¬ 
agation  constants  of  the  creeping  wave  mode  are  obtained  from  the  roots  of  the 
transcendental  equation  which  characterizes  the  coated  cylinder.  The  roots  of  the 
transcendental  equation  are  numerically  searched  via  Newton- Raphson’s  method 
by  using  both  the  Debye  and  the  Watson’s  approximations  of  the  associated  Han- 
kel  functions.  The  trajectories  of  the  roots  are  plotted  for  several  cases.  For  the 
perfectly  conducting  cylinder,  the  first  creeping  wave  mode  is  always  dominant 
over  all  other  higher  order  modes.  However  for  the  coated  cylinder,  the  dominance 
of  the  creeping  wave  mode  changes  according  to  the  change  of  the  thickness  of 
coating. 

In  the  transition  regions  adjacent  to  the  shadow  boundaries  where  the  pure  ray- 
optical  GTD  solution  fails,  a  ‘universal’  transition  integral  of  tht  uniform  GTD 
(UTD)  is  employed  to  obtain  the  uniform  solution.  Numerical  data  for  the  essential 
transition  integral  is  deduced,  via  a  heuristic  approach,  from  the  exact  eigenfunc¬ 
tion  solution  together  with  the  GTD  solution  for  the  coated  cylinder.  The  solution 
is  convenient  for  the  engineering  applications  due  to  its  simple  ray  format.  Nu¬ 
merical  results  obtained  from  the  UTD  solution  for  the  cylinder  with  thin  coating 
show  excellent  agreement  with  the  exact  eigenfunction  results. 

This  is  a  first  attempt  to  provide  a  uniform  GTD  solution  for  a  conducting 
surface  with  material  coating.  Due  to  the  numerical  difficulty  in  locating  the  higher 
order  roots  of  the  trancendental  equation,  the  examples  given  in  this  research  are 
confined  to  the  thin  coating  cases  where  either  the  first  or  the  second  mode  of  the 
creeping  wave  is  dominant.  At  this  stage,  it  is  not  yet  clear  whether  the  solution 
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is  applicable  to  the  thick  coating  case. 


For  the  furture  study,  the  three-dimensional  problem,  namely,  a  plane  wave 
incident  obliquely  upon  a  coated  cylinder  should  be  investigated.  It  is  well  known 
that,  unlike  the  perfectly  conducting  cylinder,  there  exists  both  TE  and  TM  modes 
in  the  scattered  fields,  even  though  the  incident  field  may  have  been  either  TE  ox 
TM,  The  electromagnetic  characteristics  of  the  creeping  waves  and  their  three- 
dimensional  ray  interpretations  on  the  coated  cylinder  should  be  examined.  Also, 
the  scattering  from  a  coated  elliptical  cylinder  should  be  investigated  for  the  ap¬ 
plication  of  the  solution  to  a  general  convex  surface  with  material  coating.  This 
research  is  important  and  of  interest  in  that  it  could  eventually  lead  to  an  under¬ 
standing  of  the  creeping  waves  on  the  coated  bodies  with  finite  extent. 

The  problem  of  the  radiation  from  the  coated  cylinder  is  the  case  when  the  source 
is  mounted  on  the  surface  of  the  cylinder  (  i.e.  p'  =  6  ),  and  the  field  point  is  out 
of  the  clyinder  surface  (  i.e.  p  >  &  ),  A  asymptotic  high  frequency  solution  for 
the  radiation  problem  can  also  be  obtained  by  the  similar  procedure  as  that  for 
the  scattering  problem  given  in  this  research. 
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APPENDIX  A 


WATSON’S  TRANSFORMATION 


A.l  Watson’s  Transformation 

For  the  cylinder  with  large  radius,  It  is  impractical  to  compute  the  field  from 
the  angular  eigenfunction  expansion  because  it  is  slowly  convergent.  This  difficulty 
in  convergence  can  be  overcome  by  using  the  Watson’s  transformation.  For  the 
determination  of  the  field  in  the  shadow  region  by  a  vertical  dipole  in  the  presence 
of  a  large  metallic  sphere,  Watson  (1918)  used  a  transformation  to  convert  the 
slowly  convergent  eigenfunction  series  into  a  residue  series  and  hence  a  contour 
integral  in  the  complex  plane.  He  then  showed  that  the  contour  could  be  deformed 
so  as  to  enclose  a  new  set  of  poles  and,  by  evaluating  the  residue  series  associated 
with  these  new  poles,  was  led  to  a  series  which  was  rapidly  convergent.  The 
procedure  for  converting  the  original  eigenfunction  series  (convenient  for  small 
6/Ao)  to  the  residue  series  appropriate  for  large  b/Xg  is  now  known  as  Watson’s 
transformation.  This  method  which  was  originally  used  for  the  sphere  will  be 
adapted  for  the  cylinder  of  our  interest. 

A. 2  Integral  representation 

First,  we  consider  an  integral  given  as 


-7)1 

Cx  sin  uir 


(A.l) 
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Im  [v] 


Figure  A.l:  The  contour  of  the  integral 

where  we  assume  (c^|  <  rr.  The  contour  C\  encloses  counterclockwise  all  the  poles 
due  to  the  zeros  of  sini/ff  (i.e.  t/=m)  as  illustrated  in  Figure  A.l.  Bv  in  Equation 
(A.l)  is  assumed  to  have  no  poles  on  real  axis.  Next,  the  integral  is  evaluated  via 
the  residue  theorem.  B 


zil  ,  ,  £  B(n.) 

2  ''<^1  mi^oo  ^(sin  i/7r]u=m 


(^2) 


From  (1.2),  we  obtain  the  relation 


exp\-jv{4>  ~  tt)] 


sin  t/ir 


dv  =  ^  B(m)exp(-jTn^)  (A. 3) 


assuming  that  has  no  poles  on  the  real  axis.  The  integral  in  Equation  (A. 3) 
can  be  simplified  further  if  B{u)  is  an  even  function  of  u  (i.e.  B{—u)  =  B{u)). 
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First,  we  decompose  the  integral  into  two  terms. 


/ 


^  B{u)  ^)1 

2  J-<x~j€  sin  vir 


dv  + 


^  r-<x>+jf  exp[-jV(<^  -  tt)} 
2  Joo->tj{  sin  vit 


dv 


{AA) 

By  changing  the  variable  in  the  second  term  of  Equation  (A.4),  we  obtain  the 
identity 


B(m)exp(-jm<f>)  = -j  f  ^ .  Biv) 


m=-oo 


cos  v{4>  —  ir) 
sini/TT 


dv 


{A.b) 


where  e  is  a  positive  number,  however  small.  Above  Equation  (A. 5)  represents  the 
tranasformation  of  the  eigenfunction  series  into  the  integral  form.  In  finding  the 
GO  field  in  the  lit  region,  the  integral  representation  of  the  field  will  be  evaluated 
via  the  stationary  phase  method  as  given  in  Chapter  III. 


A.3  Residue  series  representation 


On  the  conditions  that  the  poles  of  B{v)  are  located  only  in  the  lower  halfplane 
and  the  contribution  from  the  semicircle  arc  can  be  ignored  as  its  radius  approaches 
infinity,  the  integral  can  also  be  evaluated  via  the  residue  theorem  by  enclosing  the 
contour  with  a  large  semicircle  in  the  lower  halfplane,  as  shown  in  Figure  A. 2. 


Y,  B(m)exp(->m^)  =  -j  f  B{v)  dv  (A.6) 

JCx  Sin  i/TT 


m=— oo 


Hence,  the  integral  given  in  Equation  (A.6)  is  equal  to  — 27rj  times  the  sum  of  the 
residues  at  the  poles  of  B(v)  (i.e.  v^)  because  the  cuntour  C*  encloses  the  poles 
clockwise  . 


cos  1/(0  —  it) 
sin  /TT 


“oo” 

dv==2v  Y  -S(t"n) 

n=:l 


cosVn{4>  —  tt) 
sin  Vn‘>r 


{A.7) 
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Figure  A. 2:  The  contour  C*  of  the  integral 
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Thus  we  obtain 


00  **oo” 

53  5(m)exp(-jm^)  =  2?r  53  B{un)  — -r"— — ^  (>i.8) 

m^oo  „=1  sin  i/„7r 

The  above  Equation  (A. 8)  represents  the  transformation  of  the  eigenfunction  series 
into  the  residue  series  representation.  The  notation  “oo”  in  above  Equation  denotes 
in  actual  calculation  only  a  few  terms  are  needed  to  obtain  an  accurate  result.  Thus 
the  slowly  convergent  eigenfunction  series  is  transformed  into  the  fast  convergent 
residue  series.  Residue  series  solution  will  be  used  for  the  field  representation  in 
the  shadow  region. 


APPENDIX  B 


REFLECTION  COEFICIENT 

The  asymptotic  approximation  of  the  reflection  coefficient  for  the  coated  cylin¬ 
der  is  developed  in  this  Appendix.  From  the  discussion  in  Chapter  III,  the  formal 
expression  for  the  reflection  coefficient  is  given  as 


jsinag~jC^^h(t's) 
jsin  Os 


(B.l) 


In  order  to  obtain  the  approximation  of  approximation  of  should  be 

found  first. 

B.l  TMz  case 


_  Hl^hkib)d^kkia)  -  Hi^hkib)d^\kia) 
dJ\kib)Hl^\kia)  -  Hl^\kib)Hl^\kia) 


{B.2) 


From  the  Debye  approximation  of  the  Hankel  function, 


Hi^\kia)  w  y  yfe  ^  (kiasinpa  -  0akiacos(3a  -  (^-3) 


Hjj^\kia)  «  y  — - [“i  ^^iUsin/9o  -  0akiacos0a  ~  (^-4) 
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■rckibsm/3f, 


exp  ^kibs\n0h  -  0(,kibcos0i,  -  (5.5) 


H!?\kib)  «  «*P  [“J  (kibsin^b  -  /^bfeiicos/^fc  -  (5.6) 

It  should  be  noted  that  0a  ^.nd  0i  in  above  Equations  (B.3)  ~  (B.6)  are  defined 
by  the  following  relations 


V  =  kiacos0a  =  kibcos0i, 


(5.7) 


Thus,  we  obtain 


nkiJab  sin  0a  sin  0i, 


(5.8) 


Hl^\kib)Hl^\kia) 


(ab  sin/Josmp{, 


(5.9) 


where  $(a,  6,  fcj)  is  defined  as 

'4'  =  A:i[6sin/3{,  —  asin/?a]  —  ki[0f,b  cos  0i,  —  0aacos0a]  (5.10) 

Therefore,  we  obtain 

5i^^(A:i6)5i^^(/feia)  -  5i^^(/fci6)5i^^(A;ia)  « -  .}L^  sin{^(a,fe,  A:i)} 

irk\Jab  s\n0a  sin/3(, 

(5.11) 


Hl^\kib)H^J\kia)  4-  ffl^Hkib)ffi^\kia)  ^ - _i===cos{$(a,fe,fci)} 

irkiJab  sin sin^^j, 
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(B.12) 


Thus,  it  can  be  readily  shown  that 


«  jsin^b  {Hi^\kib)Hl^\kia)  +  Hl^\kib)Hi}\kia)} 


nki^Jab  sin  sin/?[> 

Therefore,  we  obtain 


cos{^'(o,  6,  fci)} 


(B.13) 


Ct  =  -j-|^sin/35Cot{^(a,6,fci)} 


(5.14) 


Incorporating  this  result  into  Equation  (B.2)  yields 

^  jsina,  -  Ij  sin/?i,cot{'f(a,6,A:i)} 
jsinQj+ ^sin/3i,cot{^'(a,6,  fci)} 


(5.15) 


In  the  Chapter  III,  it  has  been  proved  that  at  the  stationary  point. 


1/  =  kpb  cos  a,  =  kob  sin  9^ 


(5.16) 


Thus,  /3a  8Jnd  in  Equation  (B.7)  can  be  expressed  in  terms  of  the  incident  angle 

9i. 


0a  = 


cos 


(kob  sin^t 
kia 


) 


(5.17) 
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_lfkobsm6,\ 


(B.18) 


From  the  Snell’s  law,  it  is  known  that  ko  sin  =  ki  sin  6t.  As  shown  in  Figure  B.3, 
Si  is  the  incident  angle  at  the  reflection  point  and  Sf  is  the  angle  of  transmission 
into  the  dielectric/ferrite  material.  Therefore,  we  obtain 


sin/?5  =  )Jl  -  (cos/3fc)2  =  ^  1  -  ^.^cos<9,^  =  cos  Si  (5.19) 


sin  0c  =  >/l  -  (cos^a)^  ~  ^  ^  ~ 

Consequently,  we  obtain 


(5.20) 


5, 


Zq  —  jZ\  cos  ^,/  cos  Si  tan{^(a,  6,  Irj )} 
Zq  cos  +jZ\  cos  9^/ cos  Sf  tan{'i(a,6,  Icj )} 


(S.21) 


When  the  thickness  of  the  coating  d  is  very  small  compared  to  the  radius  of  the 
cylinder  (i.e.  d  4C  6  ),  ^(a,6,^l)  can  be  simplified  as 


'i{a,b,k})  kiib  -  a)sin/?{,  =  Iridcos^t  (5.22) 


If  we  use  this  result. 


j{  Zq-  jZi  cosS,/cosSi  taj\{kidcosSt)  (5  23) 

Zo  cos -fjZi  cos  ^,/ cos  tan(lrjdcos^() 

Above  result  coincide  with  the  reflection  coefficient  of  the  grounded  dielectric /ferrite 
slab. 
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B.2  TEi  case 


C  (  \  =  ~  -^  Hl^^\kib)Hl^Ukia)  -  Hl^hkib)Hl^hkia) 
^  ^  Hi^\kib)Hi^^\kia)  -  Hl^\kib)Hl^^' {kia) 

It  can  be  shown  that 

Hi'\kib)Hl^^' ik^a)  -  Hi^\kib)Hl^^' ikia) 


-j  sin l3a  {Hl^\kib)Hl^\kia)  +  AHl^\kib)Hl^\kia)} 


-4j  sin/?a 


irki  y/ ab  sin  /?«  sin  (3(, 
It  can  also  be  shown  that 


cos{V^(a,  fc,  A:j)} 


-jsinf^a  {Hl^\kib)Hi^\kia)  -  Hi^\kib)Hl^\kia)} 


4  sin  0a 


sin{'i'(o,6,  k^ )} 


irk\  ^ab  sin  0a  sin0f, 

From  Equation  (B.25)  and  (B,26),  it  can  be  shown  that 


(5.24) 


{B.25) 


(B.26) 
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(B.27) 


jsin/?6  +  H]?\kih)H^}hkia)) 


4j  sin /3a  sin/?|, 


sin{^'(a,fc,  fci)} 


irfcj  yjab  sin  0a  sin  0f, 

Substituting  Equation  (B.25)  and  (B.27)  into  Equation  (B.24),  we  obtain 


=  i‘f^sin/3i,tan{«'(a,6,*i)}  {B.28) 


Incorporating  this  result  into  Equation  (B.l)  yield 


j  sin  a,  +  ^  sin/3ttan{’i’(a,6,  fcj)} 

Rh  =  — Y -  (5.29) 

j  sinag  —  ^  sin/3{,  tan{'^(a,i>,  fci )} 


As  for  the  TMz  case  o,  and  /3j,  are  replaced  by  6^  and  St  by  using  the  relationship 
in  Equation  (B.18)  and  (B.19).  Consequently,  we  obteun 

R  -  -  jZi  cos  gt/ cos  g,tan{^(a,l>,fei)} 

^  Zo  +  jZi  cos  0tl  ^os  &i  tan{^(o,  6,  fc] )} 

For  the  small  thickness  of  the  coating  t, 


Rh  = 


Zq  —  jZi  cos  ^t/  cos  di  tan(fci</ cos 
Zo  +  3  Z\  cos  $t  /  cos  tan( d  cos  di ) 


(J5.31) 


Above  result  also  coincides  with  the  reflection  coeflicient  of  the  grounded  dielec¬ 
tric/ferrite  slab. 

In  summary, 


Zo  ~  jZi  cos$i/  co6  0t  tan{^(a,6,  fci )} 
Zo  4-  jZicos9if  cos$t  tan{^(a,&,fci )} 
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»  —  ~  cos  0i  tan{^(a,6,  fe])} 

^  Zo  +  jZi  cos  6i/ cos  6i  tan{'i'(a,  6,  fcj )} 

where 

^  =  /iri(6sin/?fc  —  asin/3a)  ~  kob{f3f)  -  /Ja)sinOi 


where 


sin  01,  —  cos  ffi 


(B.33) 


(B.34) 


(B.35) 


(5.36) 
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